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NOSEC  CARDS,  ONE  FOR  EACH  SECTION 

4  K 

THE  NUMBER  OF  WAVE  NUMBERS 
Fv/RMAT  (110) 
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Newman  has  developed  a  linearized  theory  for  the 
motions  of  a  slender  body  of  revolution,  vxth  vertical  axis, 
which  is  floating  in  the  presence  of  regular  wavas.  In  the 
present  paper  a  series  of  experimental  investigations  were 
made  and  compared  with  Newman's  theory.  Experimental 
measurements  of  motions  were  made  in  regular  and  irregular 
long  crested  waves.  Pressiytes  at  several  locations  on  the 
models  were  also  measured  ima  compared  with  the  theory. 

The  measurements  of  motions  give  excellent  agreement  with 
theory  for  slender  body.  An  extended  formula  was  developed 
for  hea*»e  motion  for  small  slenderness  ratio  of  the  body. 

The  theoretical  prediction  for  pressure  on  the  body  also 
was  found  to  give  excellent  agreement  with  the  experimental 
measurement  except  near  the  free  surface.  Observation  of 
vortex  generation  was  made  by  electrolysis. 
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NOMENCLATURE 


A 

f 

g 

T 

H 

K 

W 

*1' 


incident  wave  amplitude 
freauencv  of  oscillations  (cos) 
acceleration  of  gravity 

moment  of  inertia  of  buov  in  nitch  about  the 
center  of  gravity 

buoy  draft 

wave  number  cdVj 

radius  of  gvration  .  Fa 
buov  mass 

added  mass  of  buov  in  heave  motion 


r 

pressure 

f 

= 

ftn  = 

&) 

Section  arsA  b^oy 

S  = 

_rn_' 

w 


c  r\  =  i .  i  > 

.  1  ) 


x,y,z 

X 

e 

i 

s 

t 


To 


cartesian  coordinate  system 

nosition  of  center  of  gravity  of  bucv 

vertical  prismatic  coefficient 

density  of  flui^ 

velocity  potential 

heave  displacement  of  buov 

surge  displacement  of  buoy 

pitch  angle  of  buov 

circular  frequency  of  oscillations 

natural  period  (sec.) 


V 


K  damping  coefficient,  for  examole,  *■  *-3 

K  non-dimensional  damning  coefficient 

natural  angular  frequency 
(kft)0  the  KH  value  for  natural  frequency 


i.  intkuuUCTivjis 


The  inherent  stability  of  slender  poles  floating  with 
the  long  axis  upright  has  long  been  recognized.  Utilization 
of  this  principle  has  ranged  from  simple  wave  poles  to  more 
complex  oil  drilling  platforms.  The  concept  has  also  been 
applied  in  the  design  of  stable,  floating  platforms,  called 
"spar  buoys,"  employed  in  oceanographic  research. 

The  construction  of  a  spar  buoy  was  proposed  by  Fisher 

(31 

and  Spiess  and  reached  fruition  with  the  delivery  of  FLIP 
(Floating  Instrument  Platform)  to  the  Scripps  Institution 
in  1962.  The  original  intent  was  to  provide  a  stable  platform 
from  which  mid-ocean  acoustical  experiments  could  be  performed. 
Of  course,  the  scope  of  potential  investigations  in  the  field 
c~  -v  eanography  which  require  a  stable  platform  is  very  broad, 
~he  utilization  of  "spar  buoys"  for  various  experiments 
should  continue  to  increase  in  the  future . 

Besides  FLIP,  there  are  three  of her  "spar  buoys"  in 
existence.  One  is  called  SPAR^  (Seagoing  Platform  for 
Acoustical  Research)  and  is  operated  by  the  Navy  in  the 
Atlantic.  Another  "spar  buoy"  is  operated  by  the  Musee  Ocea- 
r.ographique  in  the  Mediterranean.  The  most  recent  addition 
to  the  growing  fleet  of  "spar  buoys"  is  called  POP^  (Perpen¬ 
dicular  Ocean  Platform)  which  is  operated  by  General  Motors 
from  Santa  Barbara,  California. 

The  interest  generated  by  FLIP  has  led  a  few  investigators 
to  develop  theories  for  the  prediction  of  the  forces  acting 
on  the  buoy  in  the  upright  direction  and  the  motions  excited 
by  waves.  These  theories,  when  practically  applied,  have  great 
value  not  only  from  the  standpoint  of  design  of  future  spar 
buoys,  but  also  for  utilization  in  oceanographic  experiments 
when  it  is  necessary  to  remove  errors  introduced  into  collected 
data  by  buoy  motions. 
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In  1963,  Newman^,  published  his  linearized  theory  on  the 

motions  of  ,!spar  buoys."  Newman's  paper  served  as  the  stimulus 

for  the  present  investigation.  It  was  used  in  the  development 

of  theoretical  predictions  for  the  pressures  acting  on  a  buoy 

in  waves  and  for  predicting  the  motions  of  a  buoy  excited  by 

waves.  Another  method  for  the  prediction  f  "spar  buoy"  motions 

12) 

was  proposed  by  Rudnick  and  may  be  viewed  as  a  simplified 
approach  to  the  problem  that  leads  to  essentially  the  same 
predicted  values  of  the  motions  as  Newman's  theory. 

The  goals  of  the  present  investigation  are  divided  into 
two  major  categories;  theoretical  and  experimental.  The 
theoretical  portion  included  a  search  for  methods  of  predicting 
the  exciting  forces  and  motions  of  a  "spar  buoy"  in  waves. 

3rtei  this  was  completed,  a  method,  based  on  the  use  of  a 
digital  computer,  was  developed  for  applving  these  theories 
to  a  "spar  buoy"  of  arbitrary  dimensions. 

The  experimental  phase  of  the  work  was  then  to  evaluate 
the  merits  of  the  theories.  For  this  purpose  a  series  of  models 
of  various  slenderness  ratios  (draft/radius)  and  bottom  ^ 
configurations  (flat  or  conical  bottom)  was  constructed.  This 
permitted  a  careful  appraisal  of  the  limits  of  applicability 
of  the  linearized  assumptions  made  by  Newman.  It  was  found 
that  the  agreement  between  theory  and  experiment  was  excellent 
for  slenderness  ratios  greater  than  about  seventeen.  However, 
as  expected,  a  modification  of  the  theory  accounting  for  added 
mass  in  heave  was  necessary  for  slenderness  ratios  lower  than 
seventeen . 
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II.  THEORETICAL  BACKGROUND  AND  COMPUTER  PROGRAM 


1.  Newman* s  Theory 

Newman  approaches  the  problem  on  the  basis  of  the  classical, 
inviscid  motion  theory.  He  seeks  a  velocity  potential,  ft*,],  Z.t) 
which  satisfies  Ladace  1  s  equation  subject  to  the  following 
boundary  conditions : 

1.  The  kinematic  boundary  condition  on  the  body. 

2.  The  free  surface  boundary  condition. 

3.  Tne  radiation  condition. 

In  deriving  the  hydrodynamic  forces  and  moments  acting  on 
the  body,  it  is  assumed  that  the  incident  waves  and  the 
oscillations  of  the  body  are  small,  and  the  body  is  slender. 

The  analysis  with  only  first  order  terms  in  the  body's  diameter 
leads  to  undamped  resonance  oscillations  of  infinite  amplitude. 

To  analyze  motions  near  resonance,  it  is  necessary  to  introduce 
damping  forces  which  are  of  second  order  with  respect  to  the 
diameter-length  ratio. 

Adopting  the  nomenclature  of  Newman's  report the 
amplitude  of  waves,  heave,  pitch,  and  surge  will  be  described 
by  A'i.V.  $  respectively  (Figure  1  portrays  the  coordinate 
system)  Denoting : 


V  —  — 

A  eHsw 

9.n0d  =  r  (  e“C*  -  ’s*)  **  Cn-  o.i  > 
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Solutions  of  Newman’s  undamped  equations  of  motion  are: 


=  rP'ai~a-(R  <<*«* 

=  ZAf - ?1&t  ~26' - ,1  Oos».t 

"  1  2  C  +Kr*-  *£)-?<  i 


From  these  equations  it  is  seen  that  resonance  occurs  in 


heave  when 


K  = 


in  pitch  and  surge  when 


K  = 


Newman  proceeds  with  his  analysis  to  compute  a  damping  term 
and  includes  it  in  a  set  of  damped  equations  of  motion.  The 
solution  of  this  set  of  equations  and  the  solution  of  the  equations 
of  the  extension  of  Newman’s  theory  to  bodies  of  small  slenderness 
ratio  ^have  been  obtained  and  included  in  the  computer  program  for 
the  evaluation  of  the  motions. 

The  solutions  of  the  damped  equations  of  motion  are  extremely 
cumbersome  and  are  not  of  practical  importance  because  computational 
experience  has  shown  them  to  be  unnecessary  for  the  configurations 
examined  except  in  a  very  narrow  region  near  the  resonant  frequency. 


defined  as  the  ratio  of  draft  to  mean  radius  of  the  body. 
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Solutions  of  Newman’s  damped  equations  of  motion  are 


K  (\-X KH& )  ^ 

sirKkyfc+tO 

5  [(i-*KH/+[i  O-XKHS.fJT 

^  -  «)* 

(5) 

(6) 

where 

bx  =  +  2V'^  -z%f 

b  ~  C-?«G«+26h) 

J|  S  »  Cft.  ~  + 

J!L_  K3  [*?,  S.8. -  &>*<?*  -  ?>6  +  Ky^  -2$f  3 
a,  =  M  CTiQi 


* 


W’i  -  ttKVi-*KM»?  I 
1  ^CXH(hXKH)  j 


P  88 

-|  =  •bamMCT^“) 
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It  should  be  noted  that  this  theory  as  well  as  the  presen- 

(2) 

tation  of  Rudnick ' s  theory  described  later  in  this  chapter  deal 
with  bodies  of  revolution  which  move  only  in  one  plane  (i.e. ,  only 
three  degree  of  freedom;  heave,  pitch,  surge) . 


2.  Extension  of  Newman's  Theory  to  Bodies  of  Small  Slenderness 
Ratio 


The  formula  in  the  preceding  section  do  not  include  heave 
added  mass.  This  is  a  consequence  of  utilizing  the  slender  body 
theory  in  the  solution  of  the  potential  problem.  As  the  model 
deviates  from  the  slender  body  assumption  made  in  Newman's  theory, 
a  modification  is  necessary. 

In  order  to  extend  Newman's  theory  r.o  include  cases  of  small 
slenderness  ratio  added  mass  and  viscous  damping  terms  should  be 
included  in  the  equation  of  motion.  The  primary  effect  of  the 
heave  added  mass  term  is  to  shift  the  resonant  frequency,  while 
viscous  damping  tends  to  decrease  the  amplitude  of  the  response 
of  the  motion  (i.e.,  to  change  the  magnification  factor). 

Since  the  damping  is  small  in  either  case,  the  effect  on 
body  motions  in  a  realistic  seaway  will  be  more  strongly  influenced 
by  changes  in  the  natural  frequency.  Therefore,  the  viscous 
damping  was  not  taken  into  account.  As  a  reference,  the  measured 
total  damping,  which  is  tie  sum  of  the  viscous  damping  and  the 
damping  due  to  the  energy  dissipation  through  waves,  is  given  in 
Table  2  in  non-dimensional  form. 

In  order  to  incorporate  heave  added  mass  into  the  equations 
of  motion,  we  first  define  two  new  added  mass  coefficients  as 
follows: 


^  >fbM' 
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It  is  well  know  that  the  added  mass  coefficient  is  a  function 
of  not  only  the  geometry  of  the  body  bat  also  of  the  frequency 
of  the  motion.  In  order  to  Find  an  approximate  value  of  this 
added  mass  coefficient  as  a  function  of  the  slenderness  ratio, 
free  oscillation  experiments  were  made  with  several  models 
representing  a  considerable  range  of  slenderness  ratio.  The 
results  of  these  experiments  are  plotted  in  Figure  5. 

The  modified  heave  motion  including  the  above  heave  added 
mass  is 


frA  (1-rKHSO _ 

Iq-XKHt  +\j  f  fjjO-XKHS.)1]2}* 


S  Joi  (O-t  -«*' ) 


(7) 


_  4-WVk(i-XkHg«,)  > 


3.  Rudnick ' s  Theory 

Rudnick  approaches  the  problem  by  considering  the  three 
major  contributors  to  the  total  force  system  acting  on  the  buoy. 
These  are  the  forces  associated  with  the  mass  of  the  buoy,  the 
hydrodynamic  accelerations,  and  hydrodynamic  response  to  transverse 
acceler'tions.  Summation  of  the  forces  and  moments  yields  one 
vector  equation  for  forces  and  one  for  moments.  Taking  the  com¬ 
ponents  and  performing  the  requisite  integrations  leads  to  the 
equations  of  motion.  For  purposes  of  comparison  of  the  two 
methods,  Rudnick* s  solutions  are  presented  here  in  terms  of  the 
same  set  of  variables  as  used  in  Newman's  theory.  The  solutions 


are: 


(7) 
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C  =  A  f  hMli]  sm  cot 

*  *  1  l-XKH  3 


$  =  Afio  c<75 


,  .  ,  Z&\  -$&> 

r  —  A  ( - 1 — S-T— )  CoS  cot 

Resonance  in  heave  occurs  when: 


K= 


w 


in  pitch  when: 


K— 


and  there  is  no  surge  resonance. 


(8) 


(9) 


y 


4.  Computer  Program 

A  digital  computer  pr^^am  was  written  in  the  FORTRAN  IV 
language  for  evaluating  the  theoretical  predictions  of  buoy  motions. 

A  listing  is  given  in  Appendix 

The  integrations  involved  in  evaluating  the  coefficients  , 

^  »  Qo  »  an<*  6-1  may  ke  performed  easily  for  several  possible 

shapes  of  the  buoy  (i.e.,  cylinderical,  conical).  These  integrations 
may  then  be  stored  in  subprograms  in  the  form  of  an  equation  for 
the  indefinite  integral  (one  subroutine  for  each  possible  shape) . 


Any  possible  underwater  configuration  of  a  buoy  may  then 
be  broken  into  a  series  of  cylindrical  and  conical  segments. 

The  sementation  of  the  buoy  is  read  into  the  main  program  along 
with  the  other  data  on  the  density  of  the  water,  the  weight  of 
the  buoy,  and  its  dimensions.  The  main  program  then  calls  on 
the  subprograms  to  evaluate  P,  ,  "P*  ,  $L0  ,  6b  for  each  segment  of 
the  buoy  and  sums  the  results  from  each  subprogram  to  obtain 
the  overall  valuer,  for  the  buoy.  Once  these  values  are  found, 
the  evaluation  of  equations  (1)  through  (9)  for  the  motions  is 
a  simple  matter.  The  actual  operation  of  the  program  is  depicted 
in  a  flow  chart  shown  schematically  in  Figure  32. 
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III.  EXPERIMENTAL  TECHNIQUE 


The  construction  of  models  and  experimental  techniques  are 
described  in  this  chapter.  Pitch,  heave,  and  surge  pressures 
and  the  vortex  generation  around  the  model  were  observed  or 
measured.  lour  methods  were  used  in  different  experiments  to 
measure  the  model  motion.  Pressure  gages  were  introduced  to 
measure  the  pressure  on  the  model  and  electrolysis  was  introduced 
in  order  to  visualize  the  flow  around  the  model.  They  are  ex¬ 
plained  in  detail  in  each  section. 


1.  Models 


Four  models  of  different  diameters  were  made  of  aluminum 
pipe.  Each  has  two  bottom  attachments.  One  is  a  flat  disk 
and  the  other  is  a  'circular  cone  shape  with  a  fixed  height  of 
three  inches  to  fit  into  the  pipe.  Therefore  the  family  incor¬ 
porates  systematic  variation  in  transverse  dimensions  only. 

The  construction  of  models  is  illustrated  in  Figure  2  and 
their  dimensions  and  characteristics  in  Table  1.  The  measurements 
performed  in  each  case  is  also  described  in  the  table. 


2.  Multiple  Flash  Photograph  Technique 

This  method  employs  a  camera  used  for  taking  still  photo¬ 
graphs.  The  camera  was  mounted  on  a  set  of  rails  about  four 
feet  long  and  placed  parallel  to  the  towing  tank  so  that  a 
picture  could  be  taken  through  the  glass  panel  at  the  side 
of  the  tank.  A  black  background  cloth  was  placed  on  the  opposite 
side  of  the  tank,  and  the  model  was  also  painted  black  with  white 
identification  stripes.  During  the  test  the  entire  tank  area  was 
from  the  stroboscopic  light  directed  at  the  model  and  set  to  flash 
at  intervals  of  one- tenth  of  a  second. 


Y 


< 


i  '‘i'.ti,  i  \  cJmuc;;  ,1  It  ;?k  t&JM  <*-*,  W/Iri  «£  Lih  1 1 
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A  test  would  consist  of  the  following  steps: 

1.  The  model  was  positioned  and  waves  generated. 

2.  All  ligh'-  were  shut  off  and  the  stroboscopic  light  was 
switched  on. 

3.  The  shutter  of  the  camera  was  opened  and  it  was  slid 
along  the  rails. 

4.  Whan  the  camera  had  moved  past  one  panel  of  the  window 
t.ie  shutter  was  closed  (about  2-3  seconds)  . 

The  film  used  was  a  polaroid  film  which  produces  both  a 
positive  and  a  negative.  The  positive  could  be  viewed  immediately 
and  the  results  either  accepted  or  rejected  after  each  test. 

If  the  results  were  satisfactory  a  large  scale  print  was  made 
from  the  negative. 

The  motion  measurements  obtained  by  using  this  technique 
were  found  to  have  a  great  deal  of  scatter.  In  addition  the 
limited  exposure  time  permitted  only  about  one  and  one-half 
periods  of  the  motion  to  be  recorded.  It;  was,  therefore  difficult 
to  obtain  a  dependable  average  of  the  amplitude  of  motion  and, 
consequently,  the  method  was  discarded. 


3.  Motion  Transducer  Technique 

In  order  to  facilitate  the  evaluation  of  results  as  well  as 
to  record  several  cycles  of  motion,  it  was  felt  that  an  electrical 
measuring  device  with  its  output  led  to  a  chart  recorder  would 
be  desirable.  Data  gathered  by  this  technique  could  be  viewed 
immediately  and  would  provide  a  record  over  many  periods  of 
motion.  In  this  was  extra..t_-oas  motions  could  be  eliminated 
leaving  only  "steady  state"  values. 

A  "constant  thrust"  gravity  dynamometer  was  available  in 
the  laboratory  and  had  previously  been  used  for  towing  ship  models 
in  head  seas.  This  device  was  too  heavy  for  the  needs  of  the 
present  experiments,  however,  so  a  new  design  was  prepared  and 


constructed  attempting  to  minimise  both  weight  the  friction  of 
moving  parts.  Figure  4  is  a  schematic  drawing  of  the  motion 
transducer.  A  description  of  this  instrument  is  as  follows. 

A  set  of  rails  was  fixed  in  the  direction  of  propagation  of 
the  waves.  A  light  weight  subcarriage  was  mounted  on  the  rails 
and  connected  by  a  wire  and  pulley  to  a  potentiometer  so  as  to 
sense  surge  translations  of  the  model.  Four  ball  bushings  were 
mounted  in  the  subcarriage  holding  two  vertical  guide  rods  which, 
in  turn,  are  attached  to  the  buoy  model  through  a  pitch  pivot  at 
their  lower  end.  Wire,  a  pulley,  and  a  potentiometer  system  is 
incorporated  between  the  vertical  rods  and  the  subcarriage  to 
measure  the  heave  motion.  The  pitch  motion  is  sensed  by  a  potentio¬ 
meter  gear  driven  from  the  pitch  pivot. 

There  are  two  drawbacks  in  this  device  which  limited  its 
successful  application.  They  both  arise  from  the  fact  that  the 
models  are  long  and  slender,  thus,  the  mechanical  connection  to 
the  model  must  be  at  the  top.  When  the  model  rotates  about  its 
center  of  gravity  both  rotation  about  the  pitch  pivot  and  surge 
translation  of  this  point,  therefore  of  the  subcarriage  result. 

The  friction  in  the  subcarriage  therefore  introduces  a  great  deal 
of  external  pitch  damping  and  at  the  same  time  the  mass  of  the 
subcarriage  increases  the  effective  moment  of  inertia  of  the  system. 
A  correction  for  these  effects  is  shown  in  Figure  7. 

The  second  detrimental  effect  is  that  tangential  motion  at 
pin  due  to  pitch  motions  must  be  subtracted  from  the  recorded 
horizontal  translations  in  order  to  obtain  the  true  surge  at  the 
center  of  gravity  of  the  model.  Since  these  terms  are  often  of 
similar  magnitude  the  errors  introduced  into  the  measured  values 
of  true  surge  may  represent  a  substantial  percentage. 


4. 


Motion  Picture  Technique 


Experimentally  this  is  the  simplest  technique  employed  in 
this  work.  A  1G  MM  movie  camera  is  set  to  shoot  through  the 
glass  panel  in  the  tank  wall.  A  clock  and  vertical  reference 
line  are  placed  in  the  center  of  the  tank  so  that  they  will  be 
within  each  frame.  The  model  is  left  free  and  positioned  so 
that  it  will  not  drift  out  of  the  frame  of  the  picture  during 
the  run.  For  length  measurement  a  reference  scale  is  marked 
c..  rhe  model  thus  alleviating  parallax  error. 

The  advantage  of  this  method  is  that  the  model  remains  un¬ 
encumbered  with  any  device  capable  of  disguising  the  actual 
motions.  The  major  disadvantage  lies  in  evaluating  the  results. 
To  say  the  least  it  is  a  time  consuming  and  most  arduous  duty. 
Nevertheless,  it  remains  one  of  the  most  accurate  and  successful 
techniques  of  measuring  the  true  buoy  motions.  It  was,  therefore 
used  as  a  "standard”  for  evaluating  the  accuracy  of  other  methods 


5.  Accelerometers  and  Rate  Gyro  Technique 

A.  Instruments  and  Experimental  Procedure 

In  the  last  two  models  (i.e..  Model  No.  3.4)  two  seismic 
accelerometers  oriented  to  sense  heave  and  surge  motions 
and  one  rate  gyro -scope  were  installed.  To  facilitate 
installation  of  these  instruments,  they  were  packaged  in 
an  instrument  module.  This  module  was  then  placed  in  the 
models. 

These  instruments  are: 

Accelerometers ; 


(for  heave) 


Statham  Lab.  Model  No.  f-2-350. 
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(for  surge?  -  Statham  Lab.  Model  No.  c-1-350, 

9  V  max  ±  1G 

Rate  Gyro 

(for  pitch)  -  u.S.  TIME  Model  No.  40,  07-90023, 

40  V.  Max.  Out  at  40  Deg/sec. 

In  this  experiment  the  r^odel  was  restrained  only  to 
prevent  it  from  rotating  about  its  vertical  axis  of  revo¬ 
lution.  This  was  necessary  in  order  to  maintain  orientation 
of  the  pitch  gyro  and  surge  accelerometer.  Two  light  alu¬ 
minum  st-uts  were  installed  to  project  upwards  on  too  of 
the  model.  The  model  was  positioned  in  the  tank  so  that 
these  struts  were  restained  between  tvro  horizontal,  parallel 
stainless  steel  rods  mounted  longitudinally  on  the  carriage 
in  the  center  of  the  towing  tank.  These  restraining  rods 
were  well  lubricated  before  aach  experiment.  The  instrument 
leads  from-  the  measuring  instruments  in  the  mcdt-l  were  made 
from  light  weight  phonograph  arm  wire  and  were  loosely  hung 
on  the  carriage  to  permit  the  model  to  move  freely  in  the 
direction  of  propagation  of  the  waves. 

It  was  initially  attempted  to  double  integrate  the 
acceleration  signals  to  get  the  corresponding  displacements 
of  motion  by  using  series  operational  amplifiers  connected 
as  integrator; .  The  final  value  was  found  to  contain  large 
errors  as  a  result  of  drift  and  noise  in  the  integrators, 
so  this  method  was  given  up.  It  was  however  found  satisfactory 
to  integrate  the  angular  velocity  output  from  the  rate  gyro 
to  obtain  the  pitch  displacement  signal.  The  accelerations  - 
the  signals  of  the  accelerometers  -  and  the  angular  displace¬ 
ment  obtained  by  integrating  the  gyroscope  output  as  well  as 
the  wave  amplitude  were  recorded  on  a  strip  chart  recorder. 


Instead  of  twic  1  integrating  the  heave  and  surge 
accelerations  the  motion  was  assumed  to  be  sinusoidal. 
In  this  case,  acceleration  is  related  to  displacement 
by  the  following  relationship, 

$(t)  =  -  (Oty)  £(*)  =  _  _L_  5ft) 


In  each  experiment  the  average  of  the  frequency  of 
the  acceleration  was  obtained  through  several  measurements 
of  accelerations  recorded  on  the  chart  paper.  It  is 
obviously  not  valid  to  use  the  above  relation  when  the 
motion  is  no  longer  sinusoidal.  However,  it  is  also 
observed ,  as  we  expected  from  the  standpoint  of  the 
shallow  water  theory,  that  the  longer  the  waves  the  more 
the  deviation  from  the  sinusoidal  waves.  It  was  net 
considered  important  to  investigate  the  responses  of  the 
motion  in  long  waves  of  which  the  wave  length  is  greater 
than  fifty  feet  since  this  kind  of  wave  is  beyond  the  expected 
range  of  large  model  motion.  For  shorter  waves  the  motion 
was  found  to  be  nearly  sinusoidal  and  this  relationship 
could  be  used. 

In  order  to  obtain  the  pure  (true)  heave  from  the 
measured  heave  motion  a  correction  should  be  applied  for 
the  effect  of  pitch.  However,  since  this  correction  varies 
as  one  minus  the  cosine  squared  of  the  pitch  angle  it  may 
usually  be  neglected.  In  the  present  case  the  maximum 
value  of  the  effect  of  the  pitch  on  heave  was  only  about 
3  to  4  per  cent  of  the  measured  heave  motion.  On  the  other 
hand,  the  effect  of  the  pitch  on  the  surge  acceleration 
measured  by  the  surge  accelerometer  was  about  90  per  cent 
of  the  measured  value.  This  was  felt  to  completely  obscure 
the  true  surge  motion  and  consequently,  surge  was  not 
evaluated  by  this  technique. 
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B.  Irregular  Waves 

Irregular  waves  are  generated  in  the  tank  by  means  of 
a  hydraulically  driven  wave  maker  in  which  the  length  of 
stroke  of  the  wave  maker  paddle  is  controlled  by  a  servo¬ 
mechanism  in  response  to  a  signal  recorded  on  magnetic  tape. 
Using  a  stock  input  tape,  the  time  scale  may  be  adjusted 
by  adjusting  the  speed  of  the  tape  recorder  and  the  amplitude 
of  the  stroke  by  adjusting  the  gain  of  the  servoamplif ier . 

The  random  waves  were  sensed  by  a  resistance  type  wave 
probe  and  recorded  in  both  digital  form  on  magnetic  tape 
and  in  analogue  form  on  a  strip  chart  recorder.  The  model 
motions  were  recorded  in  similar  fashion. 

The  digital  recorder  has  the  capability  of  sampling  up 
to  sixteen  input  channels  at  a  rate  of  one  hundred  samples 
per  second.  In  our  experiment  a  total  of  five  channels  were 
used  including  one  empty  channel  to  record  angular  displace¬ 
ment  of  pitch  motion,  heave  acceleration,  surge  acceleration, 
incident  wave  aplitude,  and  the  fifth  channel  was  empty  as  a 
marker . 

Only  one  model  condition,  the  third  model  with  conical 
bottom,  was  chosen  for  experiments  in  irregular  waves.  The 
wave  spectra  used  were  chosen  such  that  their  characteristic 
frequencies  were  close  to  the  natural  frequency  of  the  model. 

The  spectra  of  the  random  waves  are  shown  in  Figure  15. 

In  regular  waves  the  effect  of  the  pitch  motion  on  the 
heave  accelerometer  was  neglected  as  mentioned  in  the  pre¬ 
ceding  subsection,  but  this  effect  can  not  be  neglected  ir. 
the  experiment  among  irregular  waves.  For  the  pitch  motion 
is  no  lc  nger  harmonic  sinusoidal  and  the  effect  of  the  pitch 
motion  on  the  heave  accelerometer  is  not  negligible.  This 


correction  could  be  done  similarly  by  the  transformation 
of  the  coordinate  system. 

Once  the  input  (wave)  and  output  (true  motion  response; 
are  obtained  the  transfer  function  of  the  system  of  the 
motion  can  be  found.  As  output,  the  true  heave  acceleration 
and  angular  displacement  of  pitch  motion  were  used,  and  as 
input  the  wave.  A  standard  spectral  analysis  computer 
program  was  used  for  this  purpose.  In  order  to  compare 
with  this  experimental  result  in  irregular  waves,  the 
extended  form  of  Newman's  theory  was  used  as  for  a 
theoretical  prediction  of  the  response  amplitude  operator. 


C.  Calibration  and  Wave  Measuring  Device 

To  calibrate  the  accelerometers  and  integrated  rate 
gyro  output,  a  static  method  of  calibration  was  applied 
to  the  instrument  module  by  displacing  it  through  a  known 
angle.  The  angulai  output  was  thus  calibrated  directly 
and  the  apparent  acceleration  sensed  by  the  accelerometers 
was  gravity  times  the  sine  or  cosine  of  the  angle  of  in¬ 
clination. 

The  waves  were  measured  by  an  electric  resistance  wave 
probe.  This  device  consists  of  a  probe  passing  through 
the  water  surface  and  a  bare  ground  wire  on  the  tank  bottom. 
The  electrical  resistance  between  the  probe  and  ground  is 
found  to  be  very  linearly  with  the  wetted  length  of  the 
probe.  A  Wheatstone  bridge  is  used  to  detect  this  resistance 
change . 
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6.  Pressure  Measurement 


To  measure  the  pressure  samples  on  the  model  two  pressure 
gauges  were  installed  on  the  model  surface  of  the  third  model 
with  the  cone  shaped  bottom  after  the  model  as  illustrated  in 
Figure  3.  One  gauge  is  located  near  the  water  surface  and  the 
other  is  located  near  the  bottom  of  the*  cylindrical  part  of  the 
model. 

The  two  pressure  gages  used  are  as  follows: 

The  upper  pressure  gage;  Statham  Model  No.  PM  222  TC 

Sensitivity  500  MV/V/psi 

The  lower  pressure  gage;  Kulite  Model  No.  CPL-125-10 

Sensitivity  1.2  MV/V/psi 

To  calibrate  the  pressure  gages,  a  known  static  water  head 
was  applied  by  moving  the  model  with  gauges  installed  up  or  down 
a  known  distance  in  the  water. 


7.  Flow  Visualisation  -  Vortex  Observation 


It  is  of  a  considerable  interest  to  investigate  the  effect 
of  viscosity  in  the  motion  of  the  spar  buoy  in  the  waves.  In 
the  theory,  the  fluid  is  assumed  inviscid,  but  the  effect  of 
viscosity  of  water  should  be  taken  into  account  when  the  motion 
is  large  or  the  slope  of  the  body  of  the  revolution  is  very  large, 
e.g.,  near  the  edge  of  the  flat  bottom.  It  was  attempted  to 
observe  visually  the  onset  vortex  generation  near  the  flat  bottom 
of  a  model  by  introducing  into  this  region  hydrogen  bubbles  generated 
by  electrolysis.  A  fine  wire  was  put  on  the  model  and  electric¬ 
ally  isolated  from  it.  Another  electrode  was  made  of  copper  plate 
and  fixed  near  the  moving  model.  A  D.C.  voltage  was  applied  to 
these  electrodes  and  the  resulting  electrical  current  was  adjusted 
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until  the  size  and  quantity  of  the  bubbles  were  suitable  for 
observation. 

Two  different  models  were  employed  to  do  this  experiment, 
one  had  a  flat  disk  bottom  and  the  other  had  a  cone  shaped 
bottom.  The  vortex  generation  was  not  observed  in  the  latter, 
but  near  the  edge  of  the  flat  disk  bottom,  a  slight  vortex 
generation  was  observed  due  to  heave  and  large  pitch  motion. 
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IV.  EXPERIMENTAL  RESULTS 


1.  Experimental  Results  in  Regular  Waves 

The  results  of  the  motion  picture  and  motion  transducer 
measurements  of  the  experiments  performed  with  the  Model  No.  2 
with  the  cone  shaped  bottom  and  slenderness  ratio  of  about  17 
are  presented  in  Figure  6,  7,  8. 

The  resu.1  ts  of  the  accelerome  ;er  and  rate  gyro  measurements 
for  Models  Nc.  3  and  4  having  slenderness  ratios  of  about  10 
and  5  for  regular  waves  are  presented  in  Figure  9-14. 

The  added  mass  coefficient  versus  slenderness  ratio  for 
heave  motion  is  given  in  Figure  5. 


> 

2.  Experiments  in  Irregular  Waves 

The  wave  spectra  are  presented  in  Figure  15.  The  results 
of  the  heave  and  pitch  measurements  in  irregular  waves  are  presented 
as  heave  acceleration  spectra  and  pitch  angular  displacement  specula. 

The  transfer  functions  for  these  motions  were  also  obtained  and 
theoretical  prediction  of  those  transfer  function  were  computed 
as  previously  noted.  All  of  these  results  are  presented  in 
Figure  17  to  29. 

3.  Pressure  Measurement 

Pressure  measurements  for  two  different  depths  and  three 
different  angular  orientations  are  given  in  Figure  30  and  31 
together  with  their  theoretical  predictions  from  Newman's  theory. 


V.  DISCUSSIONS 


1 .  Motion  Picture  and  Motion  Transducer  Measurements 

The  results  of  motion  picture  and  motion  transducer  measure¬ 
ments  of  the  experiments  performed  with  Model  No.  2  with  cone 
shaped  bottom  presented  in  Figure  6 ,  7 ,  and  8  shows  following 
features. 

Heave .  Both  techniques  of  measurements  show  excellent 
correlation  with  the  theoretical  predictions.  The  scatter  of 
points  at  higher  KH  values  is  understandable  in  light  of  the 
experimental  error  in  measuring  these  very  small  motions  and 
the  possibility  of  introducing  at  random  exciting  forces  that 
vary  from  the  regular  wave  pattern. 

Pitch.  The  motion  picture  results  show  excellent  agreement 
with  theory.  However  the  transducer  results  are  disappointing 
and  illustrate  the  damping  which  is  introduced  by  the  instrumentation. 
The  reason  that  damping  influences  the  pitch  motion  while  4t  does 
affect  the  heave  can  be  seen  if  the  system  is  closely  examined. 

As  mentioned  previously,  pitch  friction  results  from  both  a 
rotation  about  a  pivot  and  surge  translation  along  the  rails, 
while  heave  friction  arises  from  the  tour  ball  bushings  and  the 
pulley  potentiometer  connection. 

Surge .  Again  both  techniques  yield  similar  results.  It 
should  be  noted  that  a  correction  should  be  applied  to  the  poten¬ 
tiometer  results  for  surge  in  the  case  of  Newman's  theory.  The 
effect  is  most  pronounced  in  the  area  of  predicted  surge  resonance. 

It  reduces  the  value  KH  for  resonance  slightly  and  increases  the 
value  at  higher  KH  vulues.  The  correction  is  generally  small 
and  for  the  sake  of  clarity  was  not  plotted.  The  data  does 
indicate  the  possibility  of  a  resonance  in  surge,  although  much 
more  data  is  required  before  a  conclusion  on  this  point  may  be 
reached . 
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Accelerometer  and  Rate  Gyro  Measurements 


A.  In  Regular  Waves 

The  extended  formula  of  Newman's  theory  for  heave  gives 
an  excellent  prediction  for  the  model  of  small  slenderness 
ratios. 

The  measured  pitch  motion  and  its  theoretical  prediction 
show  good  agreement  with  each  other  for  the  models  of  slen¬ 
derness  ratios  of  10  and  5. 


B .  Irregular  Waves 

The  transfer  functions  of  the  three  different  experi¬ 
mental  results  for  the  same  model  condition  show  excellent 
agreement  each  other. 

The  theoretical  prediction  and  experimental  values  of 
the  transfer  function  for  the  heave  acceleration  and  pitch 
angular  displac'anent  does  not  give  good  agreement  Quan¬ 
titatively  but  t^ey  give  good  agreement  qualitatively. 


3.  Pressure  Measurement 

The  pressure  measurements  at  the  lower  pressure  gage  near 
the  bottom  gives  a  very  good  agreement  with  the  calculation  of 
the  Newman's  theory  in  all  the  three  different  orientations. 

The  upper  gauge,  however,  shows  much  deviation  from  the  prediction 
in  all  the  three  orientations.  It  is  suggested  that  this  dis¬ 
crepancy  results  mainly  from  the  effect  of  the  free  surface, 
i.e.  we  lineexized  this  free  surface  boundary  condition. 


4 .  Applicability  of  Newman's  Theory 


For  the  model  of  which  the  slenderness  ratio  is  about  17 
Newman's  theory  gives  very  good  prediction  for  the  heave,  pitch, 
and  surge.  For  the  slenderness  of  10  and  5,  Newman's  theory 
gives  a  good  prediction  for  pitch,  but  for  heave  the  extended 
formula  of  Newman ' s  theory  gives  better  agreement  with  che 
experimental  results. 
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FIGURE  6.  Heave  Response  for  Model  #2  with  Conical  Bottom 
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FIGURE  7.  Pitch  Response  for  Model  #2  with  Conical  Bottom 
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EIGURE  8.  Surge  Response  for  Model  #2  with  Conical  Bottom 
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FIGURE  9.  Keave  Response  for  Model  #3 

with  Conical  Bottom 


Pitch  Amplitude/W.-»ve  Slope  =  f'f'l  /KA 


Newman's  Th?ory 


V*. !Situ4irsi-A»  ' 
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Frequency 

FIGURE  28.  Pitch  Transfer  Functions 
by  Newma.;‘s  Theory 


Frequency 

FIGURE  29.  Heave  Acceleration  Transfer  Fv»nc-fcTo-n.s 
by  Newman's  Theory 


FIGURE  30.  j?  vs.  KR  for  Case 


INPUT 


I  =  information  needed  tor  computations  of  subprogram 
R  =  results  of  subprogram  computations 


Possible  shapes  for  which  calculations  can  b°  made: 

(1)  Cylindrical  segments 

(2)  Conical  Segments 

(3)  Hemispherical  ends 

(4)  Elliptical  ends 

(5)  Third  degree  polynomial  end 

*  The  numbers  in  parentheses  in  the  diagram  indicate  the  steps 
involved. 


FIGURE  32.  Schematic  Flow  Diagram 
of  Computer  Program 
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APPENDIX  A:  SPAK  COMPUTER  PROGRAM  AND  ITS  WkItb-UP 


PROGRAM  SPaR  I NPiiT. OUTPUT  ) 

DIMENSION  NOSFCT ( 8 )  *  NTYPE(8)»  AZ l 9 ) »AAA( 8 ) »B8B ( 8 ) »CCC < 8 > *DDDi 8 ) 

DIMENSION  ZET  A  <  50 )  *  UNZET  A  (  50 )  »  X I  (  50 )  »  UNXI  (  50  )  ♦  PSK50) 

DIMENSION  Z£Ti  l(50) »ZETA2<  50} »ZETA3(50) »ZETA4(50) 

DIMENSION  UNPSK50) 

DIMENSION  00(50).  Ql(50).  000(50).  011(50) 

DIMENSION  AK ( 50 ) 

DIMENSION  DAMP ( 50 ) 

DIMENSION  DAMPK50) 

DIMENSION  HH( 8) 

DIMENSION  AKH{50) 

DIMENSION  PS  I  •»  (  50 )  .  XI3(50) 

DIMENSION  TITLE! 10) 

DIMENSION  EPS  I LO( 50 ) 

DIMENSION  AMOM ( 8 ) 

DIMENSION  CFREQ(50).  RPSK50) 

COMMON/FIRST/ a .  S*  C*  D*  I*  L.  SURVOL*  AREA*  TOTVOL*  VOL*  H*  CHI* 

1  W*  G»  ROE.  SO*  AZ.  AMOM.  HH*  ADD 

COMMON/SCCOND/BB*  Zl*  Z 0*  PH*  Pl»  ZG»  P22*  P2*  N*  K*  AW*  AB*  AA » 

1  CA.  Z.  A<»  000.  QO.  Oil*  01 

COMMON/ 1  HI RD/OQ*  QR •  OMEGA*  AR »  ZETA*  PSI .  XI*  EPSILO* 

1  uNZEfA*  HNPSI*  UNXI.  ZETAi*  ZeTA2*  Z£Ta3»  ZET.A4*  XI3*  PSI3 

COMMON/FOURTH/NTYPE*  AAA*  BBS-  CCC*  DDD 
COMMON/ FIFTH/ TITLE*  H£VKH»  HEVIC*  HEF*  BG»  SOAR* 

1  mover*  cfreq*  rpsi*  AKH 

C  FIRST  DATA  CARD  l‘  TUXE  CARD*  ALL  80  SPACES  ARE  READ 
READ  9.  TITLE 
9  FORMAT  (10A8) 

READ  10.  NCSEC.  G*  ROE.  W*  ZGB 
10  FORMAT  l 110.  4F10.4) 

C  NwSEC  *  NUMBER  OF  SECTIONS 
C  G  ■  ACvELERAl  IC*N  OF  GPAVITY  (FT/SFC**2> 

C  ROE  *  DENSITY  OF  FLUID  C  LB“SEc:*v?/FT**4) 

C  W  «  WPIGHT  OF  MODFL  <|  3) 

C  ZGB  «  POSITION  OF  CENTER  OF  GRAVITY  MEASURED  FROM  BASE  (FT) 

VOL  =  W/(ROE*G) 

C  VOL  *  DISPLACE!  VOLUME  REQUIRED  BY  ARCHIMEDES  (FT**3) 

C  LIMIT  PROGRAM  10  8  SErTlONS 
DO  25  I  =  1*  NOSEr 

READ  24*  NOSECT(I).  NTYPEU)*  AZ(  1+1)  »AAA(  I  J  *BBB(  1 1  *CCC  ( I )  .DDD*  I  ) 

24  FORMAT  (2110*  5F10.4) 

AMCM( I )  =0.0 

25  CONTINUE 

NOSECT ( I )  «  NUMBER  OF  EACH  SECTION  BEGlNNlNo  AT  THfc  BASE 
NTYPE  «  A  NUMBER  CLASSIFYING  THE  TYPE  0^  SECTION 

TYPE  -t  =  RIGHT  CIRCULAR  CYLINDER  R(Z»  •  CONSTANT 
TYPE  *>  =  rONICAL  SEGMFNT  R(Z)  =  LINEARLY  VARYING  FUNCTION 
TYPF  3  =  HFMI-SPHERICAL  FND 
TYPF  4  =  ELLIPTICAL  FNH 

TYPf  5  =  ARBITRARY  3RD  ORQFR  FND  R(Z»  •  A*Z**3  ♦  **Z**2  C*Z  + 
AZ ( I+t )  =  DISTANCE  FROM  THE  BASE  TO  THE  END  OF  SECTION  I 
AZ ( 1 )  =  0.0 
TOTVOL  =  0.0 
H  s  0.0 
I  »  1 


r>  r>  n  nnn  n  n  n  n  n 


5  3 


30  JOE  =  NTYPF { ! * 

A  a  AAA( I ) 

R  *  BB8<  I ) 

C  =  CCC(  !  ) 

0  a  DDD<  !  ) 

GO  TO  <720.  7->l  »  722.  723.  724).  JOE 

720  CALL  CYLINO 

c  subroutine  cylind  yields  VOLUME*  MOMENT  of  volume,  draft  and  vertical 
C  PRISMATIC  COEFFICIENT  OF  A  CIRCULAR  CYLINDER 

C  A  a  RADIUS.  B  »  C  =  D  =  0 
GO  TC  725 

721  CALL  CONE 

SUBROUTINE  CO.NL  yields  volume,  moment  of  volume,  draft  AND  VERTICAL 
PRISMATIC  COEFFICIENT  OF  A  RIGHT  CIRCULAR  CONE 
RADIUS  «  A* ( ZB"AZ ( I )  )  +  B 
A  «  SLOPE  OF  RADIUS  Vs.  2  LINE 
B  =  RADIUS  AT  AZ <  I ) .  C  =  0  =  o 
GO  TO  725 

722  CALL  HEMI 

SUBROUTINE  HEMI  YIELDS  VOLUME.  MOMENT  OF  VOLUME*  DRAFT  AND  VERTICAL 
PRISMATIC  COEFFICIENT  OF  A  HEMi-SPHERICAL  END 

a  *  radius,  b  *  c  =  =  o 

GO  TO  725 

723  CALL  ELL  IP 

SUBROUTINE  ELLIP  YIELDS  VOLUME*  MOMENT  Of  VOLUME*  draft  and  VERTICAL 
PRISMATIC  COEFFICIENT  OF  A  1/2  FLLIPTICAL  END 
A  »  MAJOR  AXIS  -  VERTICAL*  P  =  MINOR  AXIS  -  HORIZONTAL*  C  a  D  «  C 
GO  TO  725 

724  CALL  THRDOR 

C  SUBROUTINE  THRDOR  YieLOS  VOLUME*  MOMENT  OF  VOLUME*  DRAFT  AND  VERTICAL 
C  PRISMATIC  COEFFICIENT  OF  A  THIRD  ORDER  END 

C  SHAPE  Or  END*  R  a  A*Z**3  ♦  8*2 **2  ♦  C*2  +  C>  WHERE  D»0 

725  IOL  -  1 )  3C »  726.  300 

726  CALL  EXIT 

300  RFAD  301,  K 

301  FORMAT  (110) 

C  1C  *  THE  NUMBER  OF  DIFFERENT  KH'S  TO  BE  READ 
NOSEC  *  I 

READ  302.  <A<^f.)j  I  *  1.0 

302  FORMAT  (8FI0.*) 

C  A)CH  =  ICH  =  NON  -  DIMENSIONAL  NUMBER  K*H 
DO  13  I  =  :  *  < 

AMI)  a  t  <H(  I  1/H 

13  CFREO(I)  =  (G*AK< I ))**0. 5/(2.0*3.14159) 

C  AMI)  =  AW  *  1C  a  0MEGA**2/G  *  WAVE  NUMBER 

727  CONTINUE 
PRINT  7.  TITLE 

7  FORMAT  { 1H1 .  NEWMAN  SECrlCN  CALCULATIONS  *  10A8// > 

ZG  =  2GB  -  H 
Pi  a  0.0 
P?  a  0.0 
DO  303  I  =  1,  K 
OOU)  =  O.o 

303  OlCI)  =  0.0 

304  I  a  1 


n  r>  n  n  r*  n  n  n  n  n  n  n  n  n  n  n  n  no 


■j-uiu.jp 


F  4 


3U5  IF  (!  -  NOSEC >  306*  306.  360 

306  JOHN  =  NTYPPdi 
A  *  AAAI ! } 

3  =  RB8IX) 

C  *  CCCI I) 

0  «  ODDI  I ) 

IF  (AZCI+1)  -  HI  a09.  309.  307 

30?  PRINT  308 

308  FORMAT  <42H  TOUGH  LUCK  YOUR  AZd  +  T)  IS  GREATER  THAN  HI 
CALL  EXIT 

309  GO  TO  (730.  73I.  73?.  733.  734}  »  JOHN 

730  CALL  PQCYL 

SUBROUTINE  PQCtL  YIELDS  Pi*  P2»  00(N)»  Ql(N)»  (N  =  1*  K)  CIRCULAR 

SECTION 
GO  TO  30E 

731  CALL  POCC  >1 

SUBROUTINE  PC^f'N  YIELDS  Pi.  ?2»  OOIN).  Ql(N),  (N  «  1.  K)  FOR  A  RIGHT 

CIRCULAR  CONICAL  SECTION 
GO  TO  30S 

732  CALL  PQHSPH 

SUBROUTINE  PQHI  PH  YIELDS  Pi.  P2»  QOIN).  Q1(N)»  IN  «  1.  K)  FOR  A 
HEM I -SPHERICAL  END 
60  TO  305 

733  CALL  PQELL 

SUBROUTINE  PQELL  YIELDS  Pi*  P2 *  QOIN).  Ol(N).  IN  *  1»  K)  PQR  A  1/2 
ELL IPT IC *L  FNO 
GO  TO  305 

734  CALL  PQTHRD 

SUBROUTINE  PQTHRD  YIEiDS  Pi*  P2»  OOlNJ*  Ol(N).  (N  *  1*  K)  FOR  A  THIRD 
ORDER  END 
GO  TO  3C5 

360  CONTINUE 
RFAD  361.  AR 

361  FORMAT  (F10.4I 

AR  =  PAD I  US  OF  GYRATION  *  (I/M)**0.5  (FT)  ACTUAL  MASS 

READ  369,  MOVF.  "OVER 

369  FORMAT  (2110) 

IF  MOVE  =•  1  SOLVE  DAMPED  EQUATIONS 

IF  MOVE  *  2  SOLVE  UNDAMPED  EQUATIONS  (DAMPING  IK  HEAVE  ONLY) 

IF  MOVE  =  3  SOLVE  BOTH  DAMPED  AND  UNDAMPED  EQUATIONS 

IF  MCVFR  -  0  ONLY  PfPfORM  NEWMAN  CALCULAT IONS 

IF  MOVFR  *  1  Only  PEPFORM  RUDNICK  CALCULATIONS 

IF  MOVER  *  2  (OR  GREATER)  PERFORM  NEWMAN  AND  RUDNICK  CALCULATIONS 
IF  (MOVER  -  1)  370,  32*  370 

370  CONTINUE 

ADD  IS  TKE  RA-IO  BETWEEN  REAL  MASS  AND  TOTAL  MASS 
ADD  =  REAL  MA.cS  /(REAL  MASS  +  ADDED  MASS  ) 

READ  362.  ADD 
FORMAT(FlO.A) 

DO  400  I  =  1.  < 

AV  »  AK ( I ) 

OMEGA  *  ( AW*C‘)**0.5 
00  *  QO (? ) 

QR  *  Ql< ] ) 


362 


no  on 


55 


C 


IF  ( MOVF  -  2)  740*  741.  742 

740  CALL  EQDAMP 

SURROUTINF  fQo/sMP  SOLVES  THF  haMPfo  fQl'AT  1 0N5  OF  MOTION  (48)*  (49)  »  (50)  ON 

page  16  of  n^wman.  *thf  motions  of  a  spar  buoy  in  regular  waves* 

GO  TO  400 

741  CALL  FQuNOM 

SUBROUTINE  FOUNOM  SOLVES  Th-  UNDAMPED  EQUATIONS  OF  MOTION  (34).  (35)*  (36) 

ON  PAGE  i 3  OF  NEWMAN.  • ThE  MOTIONS  OF  A  SPAR  BUOY  IN  REGULAR  WAVES* 
GO  TO  400 

742  CALL  FQDAMP 
CALL  EQUNOM 

40C  CONTINUE 
32  CONTINUE 

The  FOLLOWING  IS  ATTEMPT  to  ARRANGE  A  NEAT  PRINT  OUT  OF  INFORMATION 
PRINT  4 1 C  »  TItLE 

410  FORMAT  (J.H1*  ->7H  MODEL  CHARACTERISTICS  *  10A8///) 

PRINT  31.  NOSEC 

31  FORMAT  ( 2 3H  N"MREP  OF  SECTIONS  *  110///) 

PRINT  411,  H 

411  FORMAT  (1  H  D(:AFT  =  F10.4,  6H  FT///) 

WTON  =  W/2240.0 

PRINT  33.  W,  WTON 

3 30 EORMa  r  (11H  WEIC-,JT  =  Fl2.4>  6H  LB5*3H  =  FI?. 4.  4X9HLONG  TONS// 

1  /» 

PRINT  34.  VOL 

34  FORMAT  (21)  DISPLACFO  VOLUME  =  F20.4*  9H  FT**3///) 

PRINT  40?,  CHI 

4O9  FORMAT  (4OH  VERT  I  rAL  PRISMATIC  COEFFICIENT  CHI  *  Fi0.4///) 

PRINT  35.  ZGR 

350FORMAT  ( 43H  HEIGHT  OF  CENTER  OF  GRAVITY  ABOVE  BASE  *  FlO.4* 

1  6H  FT///) 

TOTMOM  =  0.0 
DO  38  I  =  1*  -lOSEr 


38  TOTMOM  a  t{  TMOM  +  aMOMU  ) 

BCR  =  TOTMOM/VOL 
PRINT  36.  PCB 

360FORMAT  ( 44H  HEIGHT  OF  CENTER  OF  BUOYANCY  ABOVE  BASE  »  FlO.4* 

1  6H  FT///) 

PRINT  416.  AR 

4160F0RMAT  ( 37H  RADIUS  OF  GYRATION  (ACTUAL  MASS)  *  FJ2.8*  6H  FT 
1  ///) 

PRINT  418  ♦  A- 9 

4 180 FORMAT  ( 40H  RaOJUc  OF  GYRATION  (DISPLACED  MASS)  =  F\2, 8» 

1  6H  FT///) 

PRINT  200,  AR.\ 

200oformat  (ash  radius  of  gyration  about  c  (actual  mass;  *  fi2*a* 

1  6H  FT///) 

PRINT  3 ■*»  ROF 

37  FORMAT  (21H  DENSITY  OF  FLUID  =  FJ0.4*  19H  L3-SEC**2/FT#*4/// ) 

PRINT  39.  G 

3?  FORMAT  ( 2 3H  ACCELERATION  OF  GRAVITY  *  Fi0.4r  13H  FT/SEC**2/// > 

jr  (MOVFR-1)  40.  800, 4C 

40  CONTINUE 

PRINT  41,  TITL‘D 

41  FORMAT  (1H1»  13H  NEWMAN  RESONANT  FREQUENCIES  *  \0k%/ / i / / ) 


1  jjun 
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PRINT  42 

42  FORMAT  (16H  HFAVE  RESONANCE//) 

HFVKH  «  1.0/CHI 
HEVK  »  HEVKH/H 

HEF  =  (HEVK*G J**0.5/( 2.0*3.14^59) 

PRINT  43  *  HEVUF.  HEVK.  HEF 

430F0RMAT  (  7H  KH  =  Fl0.4»  10X6H  K  =  FlO.*-.  8H  FT**»1  10X 
1  6H  F  a  FlO.rt,  5H  CPS////) 

IP  <P1)  47 .  44,  47 

44  PRINT  45 

45  FORMAT  (  40H  THERE  IS  NO  RESONANCE  IN  PITCH  OR  SURGE) 

GO  TO  417 

47  PRINT  48 

48  FORMAT  ( 26H  PITCH  AND  SURGE  RESONANCE//) 

PITKH  =  Pl*H/(P2  +  AR**2  “  0. 5*Pl<r*2 > 

PITK  =  PITKH/H 

PlTF  =  (PlTK*G)**0. 5/(2. 0*3. 14159) 

PRINT  49.  PITKH.  PIT*.  PITF 

490  FORMAT  (7H  KH  =  F10.4.  10X6H  K  =  F10.4.  8H  FT**-1  10X 
1  6H  F  =  F10.4.  5H  CPS) 

PRINT  4G1,  TITLE 

4O1  FORMAT  ( IH1.30H  NEWMAN  MODEL  P.  0  VALUES  *  10A8///) 

417  PRINT  412 >  PI,  P2 

412  FORMA)  (13X6H  Pi  =  F10.4.  11H  P2  =  F10.4///) 

PRINT  413 

413  FORMAT  (6X4H  vH  5X3H  K  3X1 1H  FREO  (CPS*  7X7K  Q0(K)  7X6H  Ol(K)///) 
DO  415  I  a  l,  k 

PRINT  41/.  AKH(I).  AK(I).  CFREQd).  OOCI),  Old) 

414  FORMAT  (Flo. 2»  2F10.4.  2F15.4/) 

415  CONTINUE 

GO  TO  ( 420.  43O.  420).  MOVE 

420  PRINT  421.  TI^LE 

421  FORMAT  ( lHl »  ?7H  NFWMaN  DAMPED  MOTIONS  *  1CA8///) 

PRINT  422 

422°F0RMaT(7X.*  KH*»7X,*  K*»2X»*  FREQUENCY  (CPS>*  IX**  MAG  ZETA1/A** 
13X»*  MAG  ZETA-/A*»3X.*  MAG  XI/A*»2X»*  PSJ/KA*,?X.*  PhuCE  LAG*///) 
DO  424  I  =  1,  K 

PRINT  423*/Km(I).  AKd).  CFREQI I )  »  ZETA1  ( I  )  ’ETA2  d  ) »  XKI). 
1PSKI).  E°S1L0(I) 

423  FORMAT  (F10.2*  F10.4,  Fll.4.  5X3F13.4.  F12.4.  F12.4/) 

424  CONTINUE 
PRINT555 

555  FORMaT(1H1.*DIFFEPENT  ARRA  .GE*ENT  OF  THE  RESULT  FOR  SPECTRAL  ANAL 
1YSIS*.///) 

PRINT  666 

666  FORMAT ( 7X,<  KH*»6X»*  K*,?X.*  FREQUENCY  (CPS)**5  X.*  ZeTA3**13X.*  X 
113*.  13X,*  F Sh*///) 

444  FORMAT (F10.2»Fl0.4»4F17.4/) 

DO  1000  I  =  .,K 

lv/oO  PRINT444 »  AKH(I).  AK ( ?  i  *  CFREQd).  ZFTA31 1 )  *  XI3U).  P5l3d) 

IP  (MOVE  -  1)  440.  440.  4's0 

430  PRINT  431.  TItLE 

431  FORMAT  (  lHl »  2.9H  NEWMAN  UNDAMPED  MOTIONS  *  10A8//) 

PRINT  432 

4320FORMAT  ( 7X2HKH  7XihK  8X17HFREQUENCy  ICPS)  4X1OHMAG  ZEtA/A  6X 


r>  n 
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1  9H  f  A'i  X  I  /  A  5X10HMAG  PSI/KA//') 

PRINT  426 

426OF0RMAT  ( 5X3H  n,00  -5  ,7h  O.Onnn  8X7h  0.0000  13X7h  1.0000  8X7H  1,0000 
1  8X  7H- 1 . OnOO / ) 

00  434  I  =  It  K 

PRINT  43  3*  AKH(  I  )  »  AKU>*  CFREQ(I)»  UN^ETa(I)*  UNXKDt  UNPSlU) 

433  FORMAT  (  F 1  C< ,  2  »  flO.4*  F15.4*  5X3F15.4/) 

434  CONTINUE 
440  CONTINUE 

500  READ  501 »  I  PRES 

501  FORMAT  111!  I 

IP  (LORES  -  D  504,  600*  600 
C  IF  LPRES  *  0  MAKE  PRESSURE  CALCULATIONS  FOR  BUOY 

C  IF  LPRES  =  1  0T‘  L-nRGEP  DO  NOT  MAKE  PRESSURE  CALCULATIONS  FOR  BUOY 
504  CALL  PRSCAL 

C  SUBROUTINE  PRSCAL  COMPUTES  PRESSURE  ON  BUOY  AT  LOCATIONS  SPECIFIED 
600  CONTINUE 

IP  (MOVER  -  l<  900*  800*  800 
800  CONTINUE 

C  STATEMENTS  ROD  -  900  tRF  FOR  RUDNICK  CALCULATIONS 
8G  *  PCP  -  ZGo 
SOAR  =  AR*«2 
850  CALL  RRESON 

C  SUBROUTINE  RRESON  YIELDS  RUDNICK  RESONANT  FREQUENCIES  FOR  THE  BUOY 
CALL  RMOT 

C  SUBROUTINE  SMGT  COMPUTES  MOTIONS  OF  THE  B'JOY  USING  RUDNICK  METHOD 
900  CONTINUE 
PRINT  845 

855  FORMAT  oOH  DAMPING  COEFFICIENT///) 

PRINT  6£-6 

8560FORMAT(7X2HKH  7XaHK  8X 17HFREQUENCY  ICPS»  8X14HDAMPING  COEFF. 
1&X17HDAMP2  COEFF( ADD).//) 

DO  860  I  =  1#  K 

OAMP(  I)=W*aK(  l)*( *K(  T  )*G)**.5*«  l.“CH)*AKHC  T)*Q0(  I  ))**?/ 

1 ( 4*ROE*G*CH I **7*H**? ) 

DAMPl ( I )=\DD*DAMP< I) 

DAMP  IS  DAMPING  COEFF.  WITHOUT  THE  ADDED  MASS 
DAMPl  IS  THE  DAMPING  COEFF.  WITH  ADDED  MASS 
PRINT  857*  AKH( I ) *  AK(I).  CFREQ(I)*  DAMPl ! J ♦ DAMPl ( I) 

857  FORMAT<F10.2*F10.4*F15.4.10XF15.8*10XF15.8/) 

860  CONTINUE 
END 


SUBROUTINE  CYt.IND 

COMMON/FIRST/a*  B*  C»  D»  I*  L*  SUBVOL »  AREA*  TOTVOL*  VOL*  H*  CHI* 
1  W.  G*  ROP.  CO,  AZ*  AMOM*  HH 

DIMENSION  AZ ( 9 ) *  aMOM(8)*  HH(8) 

c  subroutine  cylind  yields  volumf*  moment  of  volume*  DR*rT  and  vertical 

C  PRISM/TI'  COEFFICIENT  OF  A  CIRCULAR  CYLINDER 

C  A  *  RADIUS.  8  e  C  =D=0 
50  AREA  =  3.14159*A**2 

S'JBVOL  =  AF  EA*(AZ(  1  +  1)  -  AZ (  I ) ) 

TOTVOL  =  TuTVOL  +  SUBVOL 
TF  (VOL  -  TOTVOL)  60.  52*  51 


i'»  i~\  r»  r> 


51  aMOM(I)  =  <0.5*'AZ<I  4  l)  -  A2  (  I )  )  4  AZ < I ) ) *SUPVOL 
1*1  +  1 

L  =  0 
RETURN 

52  H  *  AZ ( I +1 > 

AMOM(i)  =  (0,5*«AZ(I  +  1)  -  AZ(ll)  +  AZ ( I ) ) *SUBVOL 
CHI  *  W/(G*ROE*H*aREA) 

L  *  2 
RETURN 

6C  TOTVOL  =  TOTVUL  -  SURVOL 
HH( I )  =  (VOL  -  TOtVOL) /AREA 

AMOM(I)  -  (0.5*HH<I)  +  AZ(  I  )  )*3.14159*A**?*HHU> 

H  *  AZ  Cl)  +  HH< I) 

CHI  x  W/  ( G*ROE*H*aRE  A ) 

AZ ( 1+1 )  =  H 
L  =  2 
RETURN 
2ND 


AREA*  TOTVOL*  VOL*  H»  CHI » 


DRAFT  AND  VERTICAL 
CIRCULAR  CONE 


SUBROUTINE  CONE 

COMMON/FIRoT/*.  B*  C»  D»  I*  L»  SURVOL* 

1  W*  G»  ROE*  SO,  AZ.  AMOM,  HH 

DIMENSION  /Z( 9).  aMOM‘8)*  HH(fi) 

SUBROUTINE  CCNL  yields  volume*  moment  OF  VOLUME* 

PRI54ATIC  COEFFICIENT  OF  A  RIGHT 
RADIUS  *  A*(Z3-AZ(;  )  +  P 
A  x  SLOPE  OF  RADIUS  7c.  Z  LINE 
B  *  RADIUS  AT  AZ (  I  )  *  '*  =  D  =  0 

lOOOSUBVOL  =  3.14* 59*( ( ( A**2 ?/3,0?*( AZ( 1+1 )-AZ' I ) )**3  4  A*R* 
1  (AZ( 1+1 )„AZ( I ))**2  +  (B**2)*« AZ( 1+1 )~AZ( I ) ) ) 

TOTVOL  =  T(  TVOL  +  SURVOL 
DIS  *  AZ(H1)  -  AZ(I  ) 

IF  (VOL  -  TOTVOL)  110*  102.  101 
lOlOAMOMCI)  =  L«IS* (  ( 8**2  +  2.0*R*(B 

1  ( 4, 0*( R**2  +  R*(B  +  A*D IS) 

2  +  AZ  ( I  )*~UBVOL 
1*1  +  1 
L  *  0 
RETURN 

10?  H  -*  AZU+l) 

AMOM  f I )  =  DIS*((B**2  +  2.0*R*(R 
(4.0*(R**?  +  R*(R  +  A*DIS) 


+  A*DIS)  +  3.0* ( R  +  A*D IS ) **2 ) / 
4  (B  +  A*DIS)**2 ) ) )*SUDVOL 


1 

2 


A*DIS)  4  3 . 0* C B  +  A*DIS>**2)/ 
(P  +  A*OIS )**?))  )*SURVOL 


+  AZ(I)*«-UPVOL 


SO  *  3.14159* (A*(ftZr 1+1 ) — AZ C I) )  +  R)**2 
CHI  =  W/ ( G*ROE*H*cO? 

L  =  ? 

RETURN 

11C  TOTVOL  =  TOTVOL  -  SUBVOL 

HH  (  I )  =  (  ( R/ A**3 )  +  3.0*VOL-TOTVOL)/(3,14159*A**2)  )**U. 0/3.0) 
DIS  =  HH ( I )  -  AZC:  ) 

SURVOL  =  2.14i59*(  (  ( A**?)/3.0)*DIS**3  +  A*R*DIS**?  4-  «**?*r>lS) 
AMOM(I)  r  DIS*(  (?.**?  +  2.0*P*( R  +  A*DIS)  4-  3.0*(B  4-  A*DIS)**2)/ 

1  ( 4.  '**  ( fl**2  4  R*(R  +  A*DIS)  4  (  R  +  A*DIS)**2)  )  )*SUPV0L 

2  +  AZ  (  I  )*<;UBVQL 
H  x  AZ  ( I  )  ♦  HIM  I  ) 


)  =  3,  1  41  6<?»  (  A*HH'  I  )  +  R)*#2 

HI  =  w/ (r,f.R0c*H*<:0) 

'  II  +  1  )  =  W 
»  2 
•TURN 
JO 


jbroutine  hf>:i 

1MMON/FIRST/ a  ,  r.  c,  D,  j,  L.  SURVOL*  AREA »  TOTVOL.  VOL.  H*  CHI » 
w*  G.  ROE  *  50 »  AZ  »  AMOM,  hh 
mansion  azirj.  amomis),  hhiri 

JUNE  HEMI  Yielog  VOLUME*  MOMENT  OF  VOLUME.  DRAFT  AND  VERTICAL 
PRISMATIC  COEFFICIENT  OF  A  hemi-spherical  END 
vDIUS,  B  =  C  =  D  =  0 
JBVOL  =  (2.0/3.0)*3.14159*A**3 
JTVC i  =  TOTVOL  +  SURVOL 
i  'S'--  -  TOTVOL)  151,  155,  160 
,  1 'NT  16? 

I  )RMaT  (r^h  ppogram  inadequate  To  FIND  DEPTH  of  A  HEMI-SPHERE) 

['TURN 
!  =  A 

iOM(I)  =  (5«*>/S,0) *A*SUBVOL 
)  =  3,1*1 59*h**2 
<1  =  W/(G*ROF*H*cO) 

=  2 

’TURN 

HOM(I)  =  (  6,,,/8,n  j  * A# SURVOL 
*1  +  1 
=  0 
TURN 
JD 


iRROUTINF  FLLIn 

)MMON/EIR.‘ T/a,  r.  C»  0*  I,  L,  SURVOL*  AREA,  TOTVOL*  VOL,  H,  CHI. 

W,  0.  ROF,  SO,  AZ,  AMOM,  KH 
MENSIOf  AZ<9>»  aMOM(8),  HHI8I 

:>TINE  ELcIP  yields  volume,  MOMENT  Of  volume,  draft  and  vertical 
PRISMATIC  COFFFICIENT  of  a  1/2  ELLtPTlCAL  END 
,JOR  AXIS  -  VERTICAL*  B  «  MINOR  AXIS  -  HORIZONTAL.  C  -  D  *  0 
)RV0L  =  (2, 0,3,0) *3 ,141 59*»*A**2 
|  >TV0L  =  TOTVJL  ♦  SUovOL 

<VOL  -  TOTVOL)  201*  205*  210 
INT  20? 

.rmat  i 5xf2hprcc»am  inadequate  to  find  depth  or  elliptical  shape) 

TURN 
=  A/2,  C 

OM(I)  =  » 5,0/8, 0  >  * ( A/2, 0)*SUPVOL 
=  3, 14159# (R/2»0)##2 

•I  =  W/ (G*rof*h#cO) 

=  2 

TURN 

0M( I )  =  ( 5,  n/3,0)*< A/2,0 )*SUPVOL 
•  I  +  1 


J 


n  n 
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L  =  0 
RF  TURN 
FND 

SUBROUT  INF  THRDOR 

'  COMMON /F I RoT / a  »  B»  C*  0*  I*  L*  SUBVOL*  AREA*  TOTVOL*  VOL*  H»  CHI* 

1  W*  G*  ROE*  SO*  AZ*  AMOM*  HH 

DIMENSION  tZl-))*  aMOM?8)»  HH(8) 

c  subroutine  tktdor  yields  volumf*  momf.nt  of  volumf,  draft  and  vertical 

C  PR1‘«1ATI^  COEFFICIENT  OF  a  THIRD  ORDER  END 

C  SHAPE  OF  FNr>»  R  =  A*Z**3  ♦  R*Z**2  +  C*Z  +  D  WHFRE  D*0 

250  z  =  Azn+n 

SUBVOL  =  S9*Z**3*(C**?/3.0  +  Z*(R*C  +  Z*(I?.0*A*r  +  B**2) 

1  +  Z*(A*«*/3.0  +  Z*A**2/7.f')  ))) 

TmRDM  =  3.141  _>9*^**4*<C**2/4.  +  Z*«?.0*B*C/5.0  +  Z*(<2.0*A*C  +  B**2. 
1  2 )  Z6.(  +  Z*('*.0*A*B/7.0  +  Z*A**2/8.0 ) ) ) ) 

TOTVOL  =  TLTVOL  +  SURVOL 
IF  (VOL  -  TCTV0L )  251*  255.  260 

251  PRINT  252 

252  FORMAT  l SXSIHPROGPaM  INAOFOUATF  TO  FIND  DEPTH  OF  3RD  ORDER  SHAPE) 

L  =  1 

RETURN 

255  H  =  AZ  i+1) 

SO  =  3.141  R9*  I  A*AZ  (  I  +1 )  **3  +  R*AZ(|+1)**2  ♦  C«AZ(I  +  m**2 
CHI  =  W/(G*R0E*h*c0) 

AMOM(I)  =  ThRDM 
L  =  ? 

RFTURN 

260  AMOM (I)  =  THRDM 
1  =  1  +  1 
L  =  J 
RFTURN 
END 

SUBROUTINE  PQCYL 

COMMON/FIRST/A.  B»  C.  D»  I*  L »  SUBVOL*  AREA*  TOTVOL*  VOL*  H»  CHI. 

1  W*  (  *  ROF,  SO.  AZ*  AMOM,  HH 

COMMON/S rCONJ/BB*  Zl»  ZO*  Pll*  Pi*  ZG»  P22*  P2*  N,  K,  AW*  AB.  AA, 

1  CA*  Z:  AK*  000*  QO*  Oil,  Q1 

COMMON/FIFTH/tITlE?  HEVKH.  HEVIC*  HEF,  B6»  SOAR* 

1  MOVER*  CFREQ*  PPSI,  AKH 
DIMENSION  AZ ( 9 ) *  aM0M!8>»  HH(8) 

DIMENSION  A1C(  50)  •  000(50  *  00(50),  011(505  *  Ql(50) 

DIMFNSION  rFRrC-lROls  RPSK^Gj*  AKH(50) 

DIMENSION  TIT<  F  < 1 0 ) 

subroutine  pocvl  viELns  pi*  P2*  oo(N)»  gi(N)»  (n  *  i*  K)  for  a  circular 

SECTION 

310  AREA  =  3 • l4l5?*A**2 

RB  =  ARF/ *roc*g/w 
zi  =  az( i+: )  -  h 

ZO  =  AZ(  I )  -  H 

Pll  -  pB*(0.*i*(Zl**?  -  Z0**2 )  -  ZG*  ( Z 1  -  ZO)) 

Pi  =  Pi  +  Pll 

P?2  =  RB*'tZi.*3  -  ZC+*3)/3.0  -  (^1**2  ~  Z 0**2)*ZG  -+  C Z 1  -  Z0)*Z-3**?) 
i  2 ) 
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P2  =  P2  +  P22 
! c  (  I  -  1  )  702 »  70-,,  700 
^00  PRINT  70! 

701  FOR v AT  ( 1H1  ) 

702  do, sr  311.  I 

31]  FORMAT  (2?H  C  vL  I  N"R  I  CAL  SECTION  NUMBER  IlO//) 

PRINT  316.  Pi,  PH 

316  FORMAT  (6H  PI  =  Fi0.4,  11H  PAJ  =  F10.4) 

PRINT  316,  P2,  P2? 

316  FORMAT  KH  P2  =  FlO.4,  UH  P??  =  fiq.4///} 

00  3]9  N  =  1*  < 

A1--  =  A  <  { N ) 

G~3(N5  =  (RO/,  W)*<fxP(AW*Z1)  -  FXP  (  4i.'*ZD )  ) 

QO(-.i  -  0  ( N )  +  QOO(N) 

On!'!)  =  (=3R/aW)*(  FXP(  Aft*Zi  )*(Zl  -  l.C/AW  -  ZG)  -  FXP(AW*ZC)* 

1  (ZO  -  1 • a/ AW  -  ZG)) 

GMN)  =  01  (N)  +  Q,  1  (  N ) 

°R  I  NT  314.  AK.(N)»  AW.  CFRFQ(N) 

3 14"' FORMAT  (5X6H  rn  =  F10.4,  4X5H  K  =  F10.4>  8-1  FT**-1  5X 

:  13H  FfFO’iFNCY  =  -10.4,  5H  CPS/) 

=  •  I  NT  317,  N,  Qo('n,  *:00(N) 

317  •  'R?-iAf  «5H  N  =  1 5*  13H  00(N)  =  Fl0,4>  14«  COO(N)  =  FlO.4) 

PRINT  316,  N,  01CN).  Qll'N) 

318  FORMAT  (  5H  N  =  IS*i3H  Q1CN)  =  FlO.4 '  14H  Qll(N)  *  FlO.4//) 

319  CONTINUE 

I  =  T  +  1 

RETURN 

rN0 

3OOR0UTINF  POc°N 

COMMON/ F  j  R0T/-1,  o»  c »  D»  I»  L»  SURVOL*  aR'F.\,  TOT)C„,  VOL,  H*  CHI, 
1  **«  G,  ROE,  SO,  AZ,  AMGM,  HH 

COMMON/SSrONO/BB,  Zl,  ZC,  Pll,  PI,  ZG»  P22,  P2,  N,  X,  AW,  AB*  AA» 
1  CA.  Z,  A< ,  QQO,  QO,  Oil,  01 

comhon/fifth/title,  HEVKH,  HEVK »  hef,  EG,  SOAR, 

1  MOv^R,  cfr^o,  RPSI,  AKH 
DIMENSION  AZ(9),  aMOM{ 8 ) ,  HH(8) 

DIMENSION  AK (50),  000(50),  00(50),  Oij(FO)s  31(50. 

DIMENSION  CFREGI50),  RPSI ( 50 ) ,  AKH(50. 

DIMENSION  TITLE  (lO) 

C  SURROIJTINF  POCON  YJ.-'L-'S  Pi,  P2,  00  ( N )  ,  Ql(Nj,  (t-  -  1*  fQR  A  RIGHT 
C  CIRCULAR  CONICAL  SECTION 

320  Z*  =  AZ(  I  )  -  H 
71  =  AZ( 1+1)  -  H 

AO  =  B  +  \*(H  -  AZ (  I  )  ) 

BO  =  3.14159*R0E*G/W 

HU  =  oo*( A»*^*(Zi**4  -  Z0**4)/4.0  +  <2.0*A*AP  -  A*»?*ZG)* 

1  (Z  1**3  -  Z‘J**3 )  /3.0  *  (AP**2  -  ?.0*A*AP*ZG)*(ZI**2  - 

2  Z  u**-  2  )  /  2  •  3  -  AB**2*ZG*(Z1  -  ZO)) 

°1  =  PI  ♦  Pi  1 

p->?  =  qo*(a**-*(Z1**c  -  Z0**oj/5.0  +  ?.0*CA*AP  -  A**?*Z0)*{Zl**4 
1  -  Z0**4)/'.,0  +  (AR**2  -  4.0*A*AP*ZG  +  A**?*ZG**?)* 

?  (  Z 1  **3  -  Z  -»**3)  /^.O  +  2.0*(  A*ao*Zg**2  -  Ao**2*ZG) *<Zl**2 

3  -  Zr'**?’  ,?.n  +  aq**2*zg**2*( Z 1  -  ZO)) 

°2  =  P2  +  c22 


php 
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I e ( I  -  1 )  717,712,710 

710  PRINT  711 

711  FORMAT  ( 1H1 1 

712  PRINT  321*  I 

321  FORMAT  (25 H  CONICAL  SFCTION  NUMBER  110//} 

PRINT  323,  PI,  PI’ 

323  FORMAT  <6H  PI  =  Fi0.4,  11H  PJ1  =  F10.4) 

PRINT  325,  P2,  P2? 

325  FORMAT  { 6H  P2  =  Fi0.4»  11H  P'e  2  =  F10.4///) 

DO  329  N  =  1,  K 

AW  =  AK(N) 

f  R 1  NT  3?4«  AKH(N)»  AW*  C^RpQfN) 

3240F0RMAT  (  SX6H  k'H  =  F10.4,  5X5H  K  =  P10.4*  8H  FT**-1  5X 
1  13H  FREQUENCY  F10.4*  5H  CPS/) 

0  j  0  {  N  )  =  (qn/»W)*(rXP{4>.*Zl)*<A**2*<Zl**2  “  < ?.0/ AW**? ) *< AW*Z 1  - 
1  1.0))  +  ( ?.0*A*AR/AW)*( AW*Zl  -  1.0)  +  AR**+)  - 

?  EXP<  a W*Z ->)*(  ***?*( Z  0**2  -  (  ?.0/AW**2)*(  AW*Z0  -  1.0))  + 

3  < ?.0*A*A~/AW)*( AW*ZC  -  1.0)  +  An**?)) 

O'MN)  =  0n  (N)  +  0^0  (  N ) 

Of  1  <  N  J  =  (RP/aW)*(  pXP(  aW*Z  1  )  *  (  A**  2*  (  Z  1**3  -  3.0*Zi**?/aW  +  (f..0/ 

1  AW**3)*( *W*Zi  -  1.0))  +  (2.0*4*AR  ~  A**?*ZG)*(Z1**2  ~ 

2  (2.0/AW**2)*< AW*Zl  -  1.0))  +  (A 9**2  -  ?.0*A*A»*Z6)* 

3  (AW»Zl  -  l.O)/AW  -  AP**?*ZG)  “  pXP( aW*ZC)*( A**2*(Z0**3  - 

4  3.0*Z0**2/AW  +  (6.0/AW**3  )*(AW*Z0  -  1.0))  +  (2.0*A**AB 

3  -  A**2*z«r)*{Z0**?  “  (  ?.0/AW**?)*(  a  W*Z0  -  1.0))  +  (A^**? 

6  ~  ?.  '*A*;R*ZG)*( aW*ZC  -  1.0)/aW  -  AR**?*ZG)) 

GUN)  =  OICN)  *  Oi  1  (  N ) 

PRINT  326,  N,  Or ( N ) *  000(N) 

326  FORMAT  (5H  N  -  I5»13H  Q0(N)  =  F10.4»14H  000(N)  =  F10,4) 

PRINT  327,  N,  Gl(N),  Qll(N) 

327  FORMAT  (  3H  N  -  I5»13H  GKN)  =  F10.4.14H  Qll(N)  =  F10.4//) 

329  CONTINUE 

1  =  1  +  1 

RETURN 

END 

SUBROUTINE  PQHSPH 

COMMON/E I RFT/ a,  Q*  C»  D»  I*  L*  SU0VOL*  AREA*  TOTVOL*  VOL*  H*  CHI* 
1  W,  G»  ROE,  SO.  AZ»  AMOM,  HH 

COMMCN/SEC(  ND/38*  Zl*  Z0»  Pll*  PI*  Z G*  P22*  P2*  N,  < ,  AW,  AB*  AA» 
1  CA»  Z.  AK,  000,  QO.  Gil,  01 

COMMON/FIFTh/tITlp'*  HEVKH*  HEVK*  HEF,  3G,  SOAR, 

1  MOV^R,  CERpQ,  RPAI,  A<H 

DIMENSION  AZ ( 9 ) »  aMOM(8)*  HH(3) 

DIMENSION  AK(50),  Q00(50),  00(50),  011(50),  Ql(50) 

DIMENSION  CFREO(50),  RPSI ( 50 ) ,  AKM(50) 

DIMENSION  TITLE  dO) 

c  subroutine  pqhsph  yields  pi*  P2»  oo(n)»  qiin),  <n  =  i*  k)  for  a 

C  HEMI-SPHERICaL  PND 

330  R  =  A 

Ro  =  3.1415Q*R0p*G/W 
AA  s  R  -  H  -  7G 

Pll  =  RR*(2.0*AA*P**3/3.0  -  R**4/4.0) 

?1  =  Pi  +  Pll 

P?2  =  BR*(-R**5/5.0+?.0*AA*R**4/4.0  +  ( R**2  -  AA**?)*R**3/*.0 
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1  -AA*P**4  +  ;.a**2*8**3> 

P?  *  P2  +  P22 
PRINT  331.  I 

331  FORMAT  (281  H~MI -<;Phfr  ICAL  END  NUMBER  HO//) 
C  I  MUST  ALWAYS  *•  1  IF  tHIS  IS  AN  END 


PRINT 

32?  . 

PI 

»  Pll 

323 

FORMAT 

(  £  H 

PI 

=  FI 

0.4, 

11H 

Pll 

=  FlO, 

,4) 

PRINT 

325, 

P2 

,  P2*> 

325 

FORMAT 

(  6H 

P2 

=  Fi 

0.4. 

11H 

P  22 

=  FlO, 

.4///) 

DO  339 

N  = 

1* 

K 

AW  *  AK  (  N ) 

PRINT 

332, 

AK 

UN)  * 

AW, 

CFREO(N) 

332JFORMAT 

(  5X< 

H 

.'H  = 

FlO. 

4*  5X5H  K 

= 

FlO. 4* 

8H  FT** 

i  13h  frequency  =  fio.a,  5h  cps/) 

QMO(N)  =  (PB*PXP(-AW*H)/AW)*(rXP(AW*R)*(R**2  “  2.0/AW**2)  ♦  (2.0/ 

1  AW**2»*( <\W*R  +  1.0) ) 

00 ( N )  =  00 (N)  ^  000 IK) 

O')  1  ( N )  =  (ap*EXP(-AW*H)/AW)*CPXP(AW*R)*i6.0/AW**3  -  ?.0*AA/AW**2 

1  -  R**2/A^  +  AA*R**2)  “  R**3  -  3. 0* ( R**?/aW  +  (?.0/A W)* 

2  (AW*R  -  ..0)/AW)  +  AA*(R**2  +  2.0*(AW*R  +  1.0)/AW**?) 

3  +  R**?*(,W*R  +  l.O) /AW  -  AA*R**2) 

PRINT  333.  N*  QO(N).  OOO(N) 

333  FORMAT  ( 5H  N  =  I5*13H  QO(N)  =  F10.4H4K  QOO(N)  *  F10.4) 

PRINT  334.  N»  01  ( N )  *  Q1KN) 

334  FORMAT  ( 5H  N  =  I5.13H  01, N)  =  F10.4.14H  Qll (N)  =  F10.4//) 

339  CONTINUE 

I  a  I  +  1 

RPTURN 

END 

subroutine  pqfll 

COMMON/FIRST/ a.  B*  C.  D.  I.  L»  SUBVOL*  AREA*  T0TV0L*  VOL*  H*  CHI* 

1  W.  G»  roe*  SO,  aZ*  AMOM,  HH 
COMMON/SECOND /B8.  Zl,  Z0>  Pll*  Pi*  ZG*  ?**22»  °2*  N,  K*  AW*  AB*  AA, 

1  CA*  7*  AK,  000,  QO,  Oil,  01 

COMMON/FIFTH/tITLE*  HEVKH,  HEVK,  HEF,  BG,  SOAR* 

1  MOVER  *  CrRFO*  RPSI*  AKH 
DIMENSION  AZ(«J)»  AM0M(8>*  HH(8) 

DIMENSION  A<(50),  000(50),  00(50),  Q1K50),  Ql(50) 

DIMENSION  CFREQ( 50 ) »  RPSI ( 50  )  ,  AKH(50) 

DIMENSION  TITL.P  ( f  0 ) 

C  SUBROUTINE  POELL  vIELnS  Pi »  P2»  00(N)»  Ol(N)*  (N  =  1»  K )  FOR  A  1/2 
C  elliptical  END 

34 0  BR  =  3. i4159*R0E*G/W 
AA  =  H  -  \/2* 0 

CA  =  (B // )**2 

Pll  =  3B*CA*((H**4  ~  AA)/4,0  -  ( 2. 0* ( AA  -  ZG/2.0)/3.0)*(H**3  ~  AA 

1  **3 )  +  (H**2  -  H*A  -  2.0*AA*ZG)*( (H**2  -  AAJ/2.0) 

2  -  (H**2  -  H*fi )*ZG*A/2.0) 

Pi  =  PI  +  Pll 

PRINT  341 

341  FORMAT  (24)  ELLIPTICAL  END  NUMBER  110//) 

C  I  MUST  ALWAYS  -  1  *  THIS  IS  AN  END 

P?2  =  BP*CA*(  ( AA**5  -  H**5)/5.0  +  AA*(H**4  -  AA**4)/2,0 
1  (H**2  ~  W*A)*(AA**3  -  H**3)/3.0  +  2#0*ZG* ( ( <AA**4  -  H**4)/4*0 


r>  n 


7  •:  ?.0*AA*(H**3  -  AA**3>/"*.0  +  (m**7  -  H*A!**4A**7 

3  -  H**2)  /*>.  +  Zr**?*(H**?  -  h*a » *t a/7.0* n 

P7  *  P2  +  op? 

PRINT  32?-*  PI,  PI* 

323  FORMAT  ( fjH  P*  =  F^0.4»  !1h  PH  =  F10.A) 

PRINT  325.  P2.  P 2? 

325  FORMAT  <<h  P2  =  f*0.4.  I]u  P22  =  flO.4///) 

00  349  N  -  1 ,  < 

AW  *  AK(N) 

PRINT  342.  A<-«(N).  A'\,  rrRrQ<?<) 

3420F0RMaT  (  5X6H  |  H  =  FI^.a,  5X5u  <  =  rio,4,  3H  FT**-*  FX 
1  13H  pRFQmfajcy  =  r  1 0 •  4  *  CD5/  ) 

QoOlN)  =  (  FP*TA/A  v  )  *{  CXP(  AW*I-tJ  A/2.C)  )*« -AA**2  +  <?.0/A***2)* 

1  (AW*AA  +■  1.0)*(AA  -  1.0)  -  .)**->  *  H*A)  +  ?XP(-M*4W)* 

2  (2.0«H**2  +  <7.9'AW**?)*(-.\W*H  -  1.0)*{AA  -  j.O)  -  H*  a  )  ) 

00 (  N 1  =  0-  (N)  ♦  000  (N) 

0 1  1  (  N  )  =  <B3*:a/AW)  *5  gXP(-AW*AA  )  *t  AA--»*3  +  ( 3.0/4W)*4A**3  + 

1  {  6.U  ,A.V**3  )  *<  AW*AA  +  1.0)  -  ?. '1**4  -  ZO*(  AA**2  + 

7  {  7,0/AW**-?  )*<  AV»4A  +  1.0)  ;  +  {M*«7  -  u*A  -  7.0*ZO* 

3  A  A  )  *  {  /  W*  «  A  -  1  . A  )  /  A  -•  *  -  H*A  )  )  +  -XP {-***«)* 

4  <  -M  *  ■'  3  -  3.  -*H**7/AK  -  £.  •*  {  A  W*tJ  +  1.0)/AV.*r3  +  2.0* 

5  (AA  -  ZG,# i^**y  +  2« 0* l  A-*-*,J  +  1.0!  7 a.v**2 )  —  ;m**7  — 

6  H*A  -  7.  ''*ZO*AA  )*  (  AW*H  +  1.0) /A//  -  ZG*tH**2  ~  H*A  )  )  ) 

Ol(N)  =  01 (N)  +  01 j ( N) 

PRINT  343.  N,  QS(N).  DOC!  <) 

343  FORMAT  ( 5H  N  =  I5»17h  00 ( N )  =  -10.4.  l  4H  OOP(N)  =  F1C.4) 

PRINT  344,  01  Cl),  Oil  ! N ) 

344  FORMAT  (OH  N  =  1**13^  QIC’)  =  F*0.4»14H  0*1 (N)  =  F10.4//) 

349  CONTINUE 

I  =  T  +  1 

RETURN 

END 

SURROUTIN"  POTHRn 

COMMON/ FIRST  / « ,  9»  C»  D»  I»  L»  5U9V0L*  ARE*.  TOTVOL*  VOL.  H»  CHI, 
1  W*  0,  ROE.  50,  AZ,  AMOM.  HH 

COMMON/ iECOND  /3B »  Zl»  Z0»  Pll»  Pi .  Z  G»  P?7*  P?»  N,  K .  A'**  AB»  AA. 
1  CA»  Z.  A »r»  000.  00,  wll.  Cl 

COMMON/MFTh/TITlE*  HEV<H,  H  E  V  < ,  HEF,  fiG,  SOAR, 

1  MOVER*  CTREO,  RDSi,  Akh 
DIMENSION  AZ ( 9 ) ,  aM0M(8),  HH(P) 

DIMENSION  AK17.1),  QP-7J  5r- )  ♦  0^(5n),  Oll(5P),  01(50) 

DIMENSION’  CFRfo(sQ),  RPFMsG).  A<H(5n) 

DIMFNSION  TITLr  ( 1 0 ) 

subroutine  pothrd  yieids  pi.  P7,  oO(n)»  okn),  (n  =  i,  r)  fop  a  third 

ORDI  R  EN1'* 

35u  AA  *  ZG  +  H 
Z  =  AZCI+l) 

9o  =  3.14159*R0C*G/W 

P)1  =  on*Z**-a~(Z*<  '-**7/4,0  +  Z*<  7.0*n*r/S.0  +  Z*((7.0*A*r  +  o**->) 
1  /6."  +  Z*<  7.7*a*p/7.0  +  7*  A**7/n.O))))  -  a  a*  ( 

2,  +  7*19*0/?.^  Z*(l2»r>*  **r  +  d**?)/c>:-'i  +  7*  { .**o/3  .P  + 

3'  Z*A**?/7.u) ) ) ) ) 

PI  =  PI  *  Pll 

P’2  =  RH*/*  *3X(  Z**->*  (  r**2/c  .0  +  Z*(o*r/ 3.0  +  Z*({2,0*4*/"  +  0**2) 


n  n 


1  /7.G  »  Z  *■  i  A*n/4,0  +  2»A**2/9,0|)1)  -  2.0*  A  *2*  <  C**2/4 .  0 

2  +  Z*( 2.0*3*C/5.0  +  Z*(J2.0*A*C  +  R**2)/6.P  +  Z* (2.0* 

■*  A*R/7.0  +  Z*a**?/8.0) )  )  )  +  A4**?*(C**2  +  Z*(R*C/?,P 

4  +  Z*((?.n*4»-  +  R**2)/5.<">  +  Z*(A*R/3.0  +  7*A**?/Y,P )) )) ) 

P?  =  P2  +  P22 


PRINT 

323. 

PI 

.  PH 

323 

FORMAT 

(  6H 

PI 

=  Fi 0.4.  I 1 H 

Pll 

=  F 10.4 ) 

PRINT 

325. 

P2 

.  P2“> 

325 

FORMAT 

(  6H 

P2 

=  FiO.a.  UH 

P22 

=  F 1 0.4/// ) 

Z  =  -H 

+  AZ(  I 

+  i  ) 

DO  359 

N  = 

1  . 

K 

AW  *  AK  ( N } 

Z1  =  EXP(AW*Z )*( AW*Z  -  1.0)/AW**2 

Z?  =  «Z**2*FXP( AW*Z )  -  2. 0*Z 1 ) /AW 

Z3  =  (Z**3*EXP( AW*Z)  -  3.0*Z2)/Arf 

Z4  *  (Z**4*EXP( AW*Z)  -  4.0«Z3)/AW 

Z5  *  (Z**5*EXP(AW*Z)  -  5.0*Z4)/AW 

Z6  *  (Z**5*frXPC  AW*Z)  -  6.0*Z5)/AW 

Z7  *  (Z**7*FXP( AW*Z)  -  7.0*Z6)/AW 

PRINT  35? *  AKw(N).  AW.  CFRFQ(N) 

3520FORMAT  < 5X6H  KH  =  F10.4.  5X5H  K  =  F10.4.  8H  FT**-1  5X 
1  13H  FREQUENCY  =  F10.4.  5H  CPS/) 

QDO  ( N )  =  R0*<a**Z*z6  *  ?.0*A*R*Z5  +  *2.0*A*C  +  B**2)*Z4 
1  +  2.0*R*C*Z3  +  C**2*Z2) 

QD ( N )  =  OJ(N)  ♦  000 { N } 

Qll(N)  =  R0*(/.**7*Z7  +  2«0*A*B*Z6  +  (2.0*A*C  +  p**?)*Z*  +  2.0* 

1  B*C*Z4  +  r**-»*Z2  ~  ZG*(A**2*Z6  ♦  ?.0*A*B*Z5  + 

2  (2.0*/*C  +  R**?)*Z4  +  2.P*B*C*Z3  +  C**?*Z?)> 

Q1(N)  =  Q1<N)  +  Qi 1 ( N) 

PRINT  35*.  N.  00 ( N )  *  QOOIN) 

333  FORMAT  1 5H  N  =  I5»  13H  QO(N)  =  F10.4.14H  OOO(N)  *  F10.4) 

PRINT  354.  N.  01 ( N ) *  Qll(N) 

354  FORMAT  ( 5H  N  =  IC»13H  OHM)  =  FI0.4*14H  Q1KN)  =  F10.4//) 

359  CONTINUE 
1*1  +  1 
RETURN 
END 

SUBROUTINE  EODAMP 

COMMON/FIR-T/49  B*  C.  D»  I»  Li  SUBVOL *  AREA*  TOTVOL*  VOL.  H*  CHI. 

1  W,  0.  ROf,  SO,  AZ.  AMOM,  HH.  ADD 

COMMON/SECONQ/BB*  Zl,  Z0»  Pll.  Pi.  ZG*  P22*  P2*  N.  K«  AW.  AB »  AA. 

1  CA»  Z*  AK  ♦  000.  00.  Oil.  01 

common/ third /oq .  or.  omega*  ar.  zfta.  psi*  xi.  epsilo. 

1  UNZETA,  UNPSI.  UNXI,  ZETA1 ♦  ZETA2  *  ZETA3*  ZETA4*  XI3*  PSI3 

DIMENSION  AZ  (  9  )  »  aM0M(8)»  HH(8) 

DIMENSION  A<(30).  000(50)*  00(5C)*  01K50).  01(50) 

DIMENSION  4frA(50)»  PS  I ( 50 ) »  X f ( 50  )  *  fP5ilO(5P)*  UNZpTA(50)* 

1  UNPSI ( 50 J •  UNXI (50) 

DIMENSION  Z£Tfil( 50) *ZeTa2<  50) »ZeTa3( 50) *ZfTA4(50) 

DIMENSION  PS!*(50).  XI3(50) 

SUBROUTINE  EODAMP  SOLVES  THE  DAMPED  EQUATIONS  OF  MOTION  (48)*  (49)*  (50)  ON 
PAGE  16  OF  NEWMAN*  'THE  MOTIONS  OF  A  SPaR  BUOY  IN  REGULAR  WAVES* 

:;i  =  i.o- :hi*  vw*h 
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T  *  1*C  -  fH!*AW*H*00 
S  *  0,5*W*/W**V' <G*RO  ;*0MF6A ) 

TTsW*A'.v**  i*!**?/  f ->.*G*ROE*CH!*H) 

Rl  =  S*( * ?.0*Pl*GG*uR  ♦  2.0*OR**2  +  °2*00**?  -  QQ**2*R|  /  AW  + 

1  QQ**  ?<  aR**?  ) 

R?  =  Pi**?  -  ■>.->*(  ?2  -  Pl/AW  +  AR**2) 

R3  ■  2.0*(2, 0*QR  -  P 1*0Q ) 

W\  -  5*C 2.0*QQ*GR*°1  -  aR**2*QO**2  *  Pl*QQ**?/AW  -  P?*QQ**2  -  2.0 

1  *GR**?} 

W?  *  2 .0* ( f  2  _  P 1 /AW  +  aR**2 )  -  Pi**’ 

W3  »  2.0*(C<R*Pl  -  /JQ*(AR**2  +  P'  -  Pi  /AW)  ) 

C  2ETA1  IS  THE  tQlJt ION  OF  THE  EQ.  a'ITHQJT  ADDED  MASS 

C  2pTA2tI)  Io  T  hp  -0I.PT  i  ON  OF  ThF  FQUAT [ON  *ITh  ADDED  MASS 

ZEl'Al  (  I  )  =T/(Ui  **7+TT**2)**»5 
ZFTA2«  I)  =ADD»T/:ii->**2+(ADD*TT)**?)**#5 
C  PSI  =  MAGNITUDE  OF  PIt<;h  AMPL  ITUDE/wAVE  SLOPE 

P5 1  i  I  )  =  R3/< *W*(P?**?  +  R 1**2* OMEGA* *2 ) **0*  5 ) 

C  X!  =  MAGNITUDE  OF  SURGE  A*PL I  TL'DE / *’AVF  AMPLITUDE 

X  t  {  I  )  =-WA/(W->**2  +  W1**?*0MFGA**2)**0.5 

XI  (  I  )  »  ARS.'XI  (  I  )  ) 

ZFTA3  (  I  >  =  (ZE :T*2U  )*OMFGA**2)**2 
X13'.  n  =  txi  «  !  )*0MrGA**2)**2 
psi3in  =  <psi(  d*aw)**2 
EPSILC( I ) =A« 

RFTURN 

END 

SUBROUTINE  FQ.i\,<r.M 

COMMON 'FIRST/ a.  3»  C»  D*  I»  L*  oJBVOL*  AREA.  TQTVOL.  VOL.  H»  ChI » 

1  X.  C.  ROE*  SO.  AZ »  AMOM,  HH 

COMMON/^f COND/B3*  Zl*  ZO»  Pll»  Pi.  Zg»  P22»  P2»  N.  X*  A*  »  A6»  AA* 

1  CA.  7.  AK «  000,  GO,  Oil,  01 

COMMON/THIRD 'QQ.  OR.  OMEGA.  AR »  ZFlA.  PSI.  XI,  EPSILO. 

1  UNZETA,  I'NPSI.  UNXI 

DIMENSION  AZ(")»  aM0V.(8),  HH<a) 

DIMENSION  AK( 50 )  »  Q00<50),  Q0(50)»  Q11I50).  QK50) 

DIMENSION  ZETfl(30) ,  P5I<50)»  Xl<5^)»  EP~-IL0(5C)»  UNZfTAlSO), 

1  JNPSK50),  UNXI  (50) 

c  subroutine  eounqm  solves  thf  undamped  fuuations  of  motion  (34),  (35) »  (36) 

C  ON  PACE  J3  of  NEWMAN,  • ThE  MOTIONS  OF  A  SPAR  BUOY  IN  REGULAR  WAVES* 

385  'JNZETA(I)  r  ( -.  .  0  -  ch!*QQ* AW*h )  / ( 1 .0  -  CHI*AW*H) 

C  UNZETA  =  UNDAMPED  MAGNITUDE  OF  HFAVE  AMPL  I  T'JDE/WAVF  AMPLITUDE  FORMJLA  (34) 
IJNPSKI)  =-?.«*(  Pi*QQ  -  ?.0*0R)/(  AW*{  2.0*(P?  +  AR**2  -  Pi /AW)  - 
1  Pl**?>) 

C  UNPSI  ^  UNDAVPED  MAGNITUDE  OF  PITCH  AMPLI TUDE/WAVE  SLOPE  FORMULA  (36) 

UNXI ( I )  =-2.0*(Pl*GR  -  QQ*(P2  +  A«**2  “  Pi /AW ) > / ( 2.0* ( P2  +  aR**2 
1  -  PI /AW)  -  P)**2) 

UNXI ( I )  =  ARStUNXI ( I ) ) 

C  UNXI  =  UNDAMPEI  MAGNITUDE  OF  SURGE  AMPL I T JDE/ xAVF  AMPLITUDE  FORMULA  i 35 ) 
RETURN 
END 

SUBROUTINE  PRrCAL 

DIMENSION  Z(s»,  B<R(80)*  PHE(pO),  P(bO).  THETA(5) 

G=32.14 


PT=3. 14159 
RFAD10  »  R»OtL 
10  FORMAT I 2F10,4 i 
RrADl  »L  »)  »N 

1  FORMAT (  3 1 1  u  ) 

READ2* (THETA( 7  >  *  1  =  1 »L) 

2  FORMAT (8F10. 5 *» 

R£AD2*<Z( J) .J-l.M) 

003  I = 1 ♦ L 
PRINT77 

77  FORMAT (1H1) 

PRINT54.THETAM  > 

54  FORMAT!*  TKE  * NGLE  OF  THE  ORIENTATION  *  **F7.5»*  RADIANS*/) 

004  J  =  1 »  M 
PRJNT55*  Z<J> 

553 FORMAT ( // JX»*  The  DISTANCE  FROM  the  FREE  SURFACE  *  **Fl0.5»*  FEET* 
1  ) 

005  K= 1  * N 

RKR.Xj^DEL*! FLOAT <K)  ) 

ALPA*COS(TKETsU  )  ) 

BAI=2.*BKR(K)  »ALPa 
BUNE=SCRT { i»+PA!**2) 

AM*BKR(K i*L( JJ/R 
P(K)  =  EXP(AM)*n.UNE 
PhE(K)=ATAM3a!) 

C  p  IS  NON-DIMENSIONALIZEO  DYNAMIC  PRESSURE»P/RHO*G*A 
5  CONTINUF 
PRINT  1000 

loOO  FORMAT {//15X»*NO.*»17X**KR*»12X**P/RhO,G*A*»9X»*PhASE  ANGLE*/) 
00100  K- 1  »N 

iv-C  PRINT99*  8*R<K).  P(K),  PME(K) 

99  FORMAT (8X • 1 10  .3F20.4/ ) 

4  CONTINUE 

3  CONTINUE 
RETURN 
FND 

SUBROUTINE  RRESON 

COMMON/FIRST/ a  *  B»  C»  D*  I*  L»  SUBVOL*  AREA*  TOTVOLt  VOL*  H*  CHI* 

1  W.  G*  ROE*  SO*  A Z.  AMOM*  KH 

COMMON/SECONO/BB*  Zl*  Z0»  Pll»  Pl»  ZG*  P22*  P2*  N*  K»  AW.  AB*  AA* 

1  CA*  Z *  AK*  QOO.  Q0»  Qll*  01 

COMMON/F IFTH/tITLE*  HEVKH.  HEVK.  HEF,  8G*  SOAR* 

1  MOVER.  CFRFO,  RPSIf  AKH 

DIMENSION  ',Z(9>»  &MOM(  8  )  »  HH(8) 

DIMENSION  CFRrO(50)»  RPSI ( 50 ) *  AKH(50> 

DIMENSION  AX(.«.0)  .  000150).  Q0(  50)  *  Qll(50j»  Ql(50) 

DIMENSION  TITLE  C 1 0) 

C  SUBROUTINE  RRESON  YI£ldS  RUDNICK  RESONANT  FREQUENCIES  FOR  THE  BUOY 
PRINT  1.  TITLE 

1  FORMAT  <ihi*  34H  PUDNICK  RESONANT  FREQUENCIES  *  \0^/////) 

PRINT  3 

3  FORMAT  ( 16H  HEAVE  RESTNANCE//) 

PRINT  7.  HEVKH.  HEVK.  HEF 

70FORMAT  ( 7H  KH  =  FlO.4*  10X6H  K  =  F10.4.  8H  FT**“1  10X6H  F  * 


wrssre>) 


1  FiO.4,  5H  C?S////) 

PRINT  10 

FORMAT  (16H  PITCH  RE SONANCF/ / ) 

PTTK  =  Pi/(P2  +  SOAR) 

P I TKH  =  pitk*h 

PITF  =  (  PI  TK*G)  **0.5  /  (  2.  0*3*  ■>  4159 ) 

PRINT  15*  PITKH*  P I T< *  PITF 
I50F0RMAT  (7H  <H  =  F10.4,  10X6H  <  =  FIG. A 

1  F10.4*  5h  CPS////) 

PRINT  20 

FORMAT  (  28H  ThERf  IS  NO  SURGE 

RP  TURN 

END 


<  =  F10.4*  8h  FT**-2  10X6H  F  = 


SUBROUTINE  RMOT 

COMMON/Fj«oTM.  B  ♦  C »  D»  I*  L.  GURVOL*  A«EA»  TOTVOL*  VOL*  h»  CHI* 
1  W*  <}»  ROr  •  SO*  A 2*  AMOM,  HH 

COMMON/S LCOND'BB*  l\ .  20,  Pll*  Pi*  2G*  ?22»  P2*  N,  < *  Ah •  AB*  AA, 
1  CA»  2*  A;.»  Q00 *  Q0»  Oil,  01 

common/third.-oo*  OR*  omega*  AR»  2eTa»  PSI*  Xr,  EPSTLO* 

1  UN2FTA*  "NPSI*  UNXJ 

COMMGN/FIFTH/T I TlE »  HEVKn*  H£VK»  nEF*  BG»  SOAR* 

1  MOVER*  CFRFO*  RPSI*  AKH 

DIMENSION  A2(,')»  a,M0M(8)»  HH(8) 

DIMENSION  /<(  '0)  ,  000(50),  G0(50),  01K50),  0K50) 
^IMENSl0N^(ErA(50,*iPSH50,,  X,.5.„  FPSlL0(50,,  UN2£Ta<50), 

DIMENSlOf  CFRE0(50)»  RPSI ( 50 ) »  AKH(50) 

DIMENSION  TITl-  ( T  0 ) 

6UBR0UTINE  COMPUTES  MOTIONS  OF  THE  8U0Y  USING  RUDN1CK  METHOD 

F’KINi  4*  TITLF 

1  FORMAT  (1H1*  ->  1H  PUDNICK  MOTIONS  *  10A8///1 

PRINT  3  1  AB///) 

30r°RM*T  ,7*-HKH  7X,HE  exlTHFREQOENCY  ICPS)  *X]OHMAG  ZETA/A  6X 
1  9H  MAO  XI /A  7X  10HMAG  PSI/KA///) 

PRINT  5 

5GF0RMAT  (5XSH  a.-«0  3X7H  0.0000  8X  7h  0.0000  13X7H  l.noop  8X 
1  7H  l.OOOo  8X7H-1.0000/) 

IF  (MOVER  -1)7,15,15 

7  PRINT  8 

8  callAex(?5H  EXIT  PM0T  STEP  8  M0VER  LFS"'  THAN  1» 

15  CONTINUE 

DO  25  i  =  1,  K 

=  (2,A*Ql(I)  ~  P1*00(  I  )  )/(  aK{  I  )*{P2  +  SOAR-  Pl/AKd))) 

1R  SiSIr  AKUl’  CFREQII>’  ONZETAIl),  QOU).  RPS I C i! 

18  FORMAT  (Flu. 2*  F10.4*  F15.4,  5X3F15.4/) 

25  CONTINUE 

return 

END 
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J'S~- 


WAVE  NUMBER  TIMES  DRAFT  *  KH  (FT) 

FDRMAT(8Fl-.45 

8/K  CAROS  -  8  KH.S  PER  CARD 

6  AR 

RADIUS  OF  3YRATION  OF  ACTUAL  MASS  (FT) 

FORMAT (2Flu. 4 i 
1  CARD 

7  MOVE -MOVER 

numbers  TELLING  MA:_  «!N~  ifc*HICH  COMPUTATIONS  TO  PERFORM.  (SEE  NO. 
369  ) 

FORMAT (2110) 

1  CARD 

3  ADD 

THE  ADDED  MASS  COEFFICIENT  «  THE  RATIO  OF  REAL  MASS  AND  TOTAL 
MASS 
1  CARD 

9  L^RFS 

number  tflling  MACHINF  WHETHER  to  compute  pressures 

SPE  NO.  50U 
FORMAT ( I 1 - 1 
1  CARD 

***  FOR  PRESSURE  CALCULATION  *** 

1^  R.DEL 

MAX  RADIUS  OF  MODEL (FT).  INCRFMFNT  OF  KR  VALUE 
FORMAT ( 2Fll'.4 ) 
l  CARD 

11  l*m»n 

NUMRPR  OF  orifntations.numrfr  OF  dfpths.  number  OF  KR  VALUES 

FORMAT ( 3 1 1U ) 

1  CARD 

12  THETAd)  (  I«1  »L) 

ANGLE  OF  ORIENTATIONS  (RADIANS  ) 

FORMAT 0*10. 5) 

L/8  CARDS 

13  Z(J)  (  J>»  1  *  M  ) 

DISTANCES  FROM  FREE  SURFACF  (FT) 

FORMAT ( 8FlU,5 ) 

M/8  CARDS 

*  this  is  the  last  data  card 
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r  sample  pro*lcm  input  f'ATA 


’O^fL  A  CONE  POTTT'  (?  ^  TXP.) 

2  1.9^67  11.78  .7<>1 

1  7  .78  *  ’’E 

7  1  7.  .1878 

8 

.5  l.~  1.8  2.0  7.5 

.64? 

3  3 
.P4  78 

.1876  l.l 

3  7  ]  U 

.0  1,571-3  3.141*9 

•178  — 1 .48*3 


CONIC  L  i»t  C T 1  Civ  NtiMdfcW 


t 

f 


i  I 

i  } 


in 

<0 

in 

Si 

Irt 

a 

u 

a. 

GL 

a 

u 

o 

o 

U 

u 

♦ 

(V 

01 

cr 

1 

r-  i 

n 

9! 

Ol 

o 

o 

o» 

SI 

o 

*1« 

a* 

f 

o 

9 

o 

• 

• 

• 

• 

• 

a 

a 

a 

a 

a 

>- 

>- . 

>> . 

>- . 

>-: 

o 

o 

o 

u 

u 

z 

z 

z 

z 

z 

oJ< 

oJ 

oJ 

Ui 

U 1 

a 

3 

3* 

3 

3 

3f 

31 

SI 

34 

31 

ol 

il 

xX 

'.U 

ol 

X* 

9 

>*■  • 

X 

»  > 

x- 

o. 

91 

X* 

O  s> 

X 

■*-• 

J.I 

-« 

OI 

i.t 

>«  -• 

0-1 

M* 

r»-  ■ 

J.I 

9!#< 

o-l 

i>  m 

** 

~4 

— « 

o 

o  o 

o  o 

3 

z> 

O  3 

3 

o 

o  O 

o  o 

• 

♦ 

•  • 

• 

• 

•  • 

•  • 

t 

i 

1 

* 

• 

. 

*•* 

—4 

_ 

1 

1 

1 

1 

1 

1 

f 

a 

• 

♦ 

a 

a 

♦ 

a 

M 

♦ 

a  a 

• 

a 

a 

a 

a  a 

a 

a  a 

r- 

»- 

H- 

a 

m. 

■M> 

U» 

— * 

*A-  1 

X 

-9  -T 

z 

•9 

z  z 

z  z 

— » 

w- 

■— ■ 

© 

o 

o  ^ 

o 

o 

«-4 

3 

O  — • 

o 

O  r-4 

« 

3 

— # 

O  H 

9 

3 

-*• 

♦ 

o  -*• 

o 

3  -«• 

9 

3 

3 

.n 

o  J 

> 

3 

3 

9 

3  3 

n 

3  3 

Ol 

j) 

r- 

3 

*»  9 

• 

• 

• 

• 

• 

7*  ^ 
M  -• 

**  o 


a  a 

^  j' 
-•  M 
X  X 


o  c 

r 

v 


a 


s>  i- 

-•  30 

n  *•+ 

o  o 


•X)  •> 

o 


a 

z 


o  o 

oj  r*- 


o 

i 


O  9 
»  * 


a 


a 

c 


m  •** 
o  o 
o  o 


U 

II 

II 

a 

a 

a 

a  a 

a 

a 

9 

3 

o 

9 

9 

9 

3 

9 

9 

— *-» 

3 

O 

Z 

Z 

3 

Z 

z 

O 

z 

3 

z  z 

O 

z 

z 

IT) 

C 

m 

— * 

3 

^ « w 

in 

— * 

• 

-* 

• 

9 

—4 

• 

o 

«-4 

• 

3  -• 

• 

c 

—4 

3 

3 

-4 

3 

3 

-4 

3 

3 

M 

3.  3  . 

\i 

3 

3 1 

If 

- 

n 

,\i 

a 

u 

a 

a 

a 

in 

x. 

r 

r 

X 

X 

r 

V 

5c 

5C 

a 

a 

it 

u 

a 

a 

a 

a 

it 

a 

■» 

z 

z 

z 

z 

Z 

z 

z 

z 

z 

CYl  I  vi^iCoL  '  1  u-'x  .iiror- 


•m 


^■PBWPBP 


I 

72 


</) 

IA 

(A 

(A 

Ul 

a 

a 

a. 

a 

a 

o 

o 

u 

u 

u 

*• 

M 

M 

o 

r*. 

d 

C: 

O 

o 

3 

o« 

# 

• 

n 

o 

• 

o 

■D 

m 

*>• 

3 

* 

O 

•  • 


VI 


t  t 
VI 
DOD 
MM 
•  • 


t 


IS  X 

—  r 

—  \j 

i.  x 


\J 

O 


N 


U 
2 
iJ 
O 
3 1 
J 

r- 

4.1 


I 


v/ 

o 

o 

VI 


.3 


n 


o 

o 

3 


a 

a 

K 

>- 

► : 

»-. 

o 

SJ 

(J 

z 

Z 

z 

aj 

u 

oli 

n 

3 

3 

3 1 

31 

3  ( 

ii 

il 

ul 

n 

z- 

VI 

til 

r  ■ 

n 

r  « 

D 

x> 

. 

i.  I 

o. 

O* 

i.t 

VI 

a.  1 

A 

>n 

*n 

♦ 

*» 

•n 

i) 

(M 

\i 

o 

M 

n 

M 

* 

M 

• 

m 

• 

• 

• 

• 

• 

1 

> 

1 

♦ 

* 

♦ 

»* 

* 

a 

M 

♦ 

N 

W 

♦ 

m 

n 

*— 

»— 

*— 

j* 

a. 

■» 

a. 

•* 

a. 

* 

- ' 

tl. 

7, 

O 

o 

3 

O 

-a 

o 

o 

—4 

-a 

3 

— • 

D 

3 

— « 

* 

3 

-a 

3 

3 

3 

C* 

3 

3 

O 

3 

3 

** 

• 

• 

M 


^  . 

U 

Z 

Ul: 

3 

31 

U 

X 

i.1 


I 

* 


o 

o 

••n 

■*» 


c 

i*j  O 
o  r*j 
»  D  M 

•  • 

C 


cr  > 

•o  -a 

n  n 
a  r  M 

•  • 

c 


3  VI 

*  c 

♦  ■*> 
■  O'  (V 

•  • 

£. 


S3  ~ 

■>  o 

A  VI 

M  M  H 

•  »i 

<■  Z 


U  M 

•3  3 

3 -  -3 

S3  Z  z  o 

C  —  —  ;f> 

•  • 

»l  J  3  vj 


»  n 

—  A 

-  o. 


M  M  V 
C 

v 


II  -n  ■»> 

r 


x  *  * 

r 

y 


n  ji  r 

y 


n  ii 


M  N 

Z  Z 


a  ii 


U  N 


N  U 

z  z 


z 


*?!?*  *  (mU?C1  7  0  T  i»  •  a  (M)Ta 

/<*«**'  *  (fc'MJftrt  OfftC*  «  ( M )  ^  r> 


73 


wCuEl  vnartACTfc-HlS' ICS  •  MODEL  A  CONc  aPlTUW  ip  ND  EAp.) 
NUMhEo  0?  SeCUunG  ■  2 

DRAFT  *  l«B79t>  FT 

WEIGHT  *  U.7d00  <rp  *  .0053  LONG  TONS 

DISPLACED  VOLUME  *  *1092  FT**3 

VERTICAL  PRISMATIC  COEFF ;C1EnT  *  C?I  ■  *9113 

HEIGHT  OF  CENTtH  OF  GKAV J  T*  ABOVE  BASE  a  #7900  F] 

height  of  center  of  buoyancy  above  base  «  1..0221  ft 

radius  Of  GYRATION  {ACTUAL  NASS)  *  .6*200000  FT 

radius  of  gyration  (uisplaeeq  mass)  »  I  ft 

RADIUS  of  GyRaIION  ABOUT  C  (ACTUAL  MASS)  ■  I  FT 

DENSITY  OF  FLUu  *  1*9367  L8-SEc**2/FT*#4 

ACCELENaTIO».  C£  GRAVITY  a  32.1300  FT/SEC**2 
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*OUEL  U  VALUES  *  MUutL  a  CG(4  BOTTOM  (2  NO  e*P,) 


Pi  a 

.232]  P2  * 

•  3002 

Kn 

H 

FREn  (CPS) 

Oft(«M 

Ui  (K> 

•3l> 

•  £660 

.40&4 

.6030 

•  C3O0 

1  «u0 

?b32« 

.63*2 

•  6*>6o 

•  es**i 

1  «30 

•  7980 

.B0O2 

.5449 

•  2302 

5,i»0 

1  »C64u 

.9JU9 

•  43*8 

•  2207 

?«3C 

1.33C* 

1.C4U7 

•  3^8 

•  fiol 
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rwe  AMtiLfc  OF  fn£  0**Ii«TATiyN  *  0.  «*0IANS 


rHt  ol«T>*tE  FkO*  T^£  FR££  SURFACE  ■  -.12-90  FtET 


*0, 

AR 

P/RHO.O.A 

FMASE  A*6LE 

1 

•  1000 

.9540 

•  1974 

2 

•  2000 

.9426 

•  3005 

3 

•  3000 

•  9548 

,6404 

4 

•  4000 

•  9009 

•  0747 

5 

•  5000 

1.0153 

.7054 

6 

•  6000 

1.0*71 

.0761 

7 

.7000 

1*0709 

.9505 

a 

.8^00 

1.1069 

1.0122 

q 

.9000 

1.1301 

I,0o37 

10 

1.0000 

1.1*80 

1.1071 

Trtt  Ol«:Tfc<Ce  FKv p  THt  F«EE  SURFACE  a  -1,45830  FfcfcT 


Mi. 

*R 

p/RHO.o.A 

RMASE  AN61.E 

\ 

•  1000 

•4M5 

•  1.74 

2 

•  2000 

.2*73 

•  5805 

3 

•  3000 

•  llfl 

,8404 

4 

•  4060 

•  0&?i 

.6747 

5 

•  5000 

•  0289 

,7854 

6 

•  6000 

•  0147 

•  6761 

7 

•  7000 

•  0074 

.9505 

a 

•  8000 

•  0057 

1.0122 

9 

•  9000 

•  0019 

1.0637 

10 

1,0000 

•  0009 

1.1071 

! 
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THE  AW(*LE  OF  THt  ORIENT ATiy*'1  *  l«5f080  R*UlANS 


T*t  Dls^-^CE  FkO* 

T-4  FREE  5U?FACE  a 

-.12500  FttT 

vo. 

HR 

p/RMU.O, A 

PRASE  AIV6LE 

1 

•  1000 

•  9355 

••0000 

? 

.200  0 

.8752 

•*0000 

3 

.3000 

•  8187 

-•0000 

4 

.4000 

•  7659 

-.0000 

5 

•  5000 

.7165 

-.0000 

6 

•  6000 

.6704 

— • oOQO 

7 

•  7000 

•  6f  7 1 

-.oooo 

a 

•  8i/00 

•  5566 

-.UUOO 

9 

•  9000 

•  5466 

-.oooo 

1  0 

1  «0o00 

•  5134 

-.0000 

T«c  nIsT*'vCE  £  kO* 

ThE  Free  SURFACE  a 

-i ..45530  Fttr 

HU  . 

ISR 

p/RHO *Oa A 

FnASE  aiv^LE 

1 

•  1000 

•4594 

-.oOOO 

? 

•  2000 

.2111 

-•OOOO 

3 

•  3000 

•  0970 

-.0000 

4 

•  4000 

•  0446 

..oooo 

s 

.5000 

•  0305 

— • oooo 

6 

•  6000 

•  0094 

-•oooo 

7 

.7000 

•  0044 

-.0000 

8 

•  6000 

•  0030 

-.0000 

9 

•  9uQ0 

.0009 

—.0000 

19 

1.0000 

•  0004 

-.0000 

81 


THE  Of  (*t  0*iCNTATlOH  ■  j.i*ls%  riuIams 


THE  OlsTA'tCE  7HO* 

the  free  surface  . 

-.12600  FfetT 

•<0. 

KR 

P/RHO.tt.A 

Phase  ang^e 

X 

•  1000 

•  9640 

-•1?74 

2 

•  2000 

•  9426 

— • 8605 

3 

•  3000 

.96*0 

-•5404 

4 

•  40Q0 

•  9809 

••$747 

5 

•  5000 

l?0l33 

-•?854 

6 

•  6000 

1.0471 

-.8761 

7 

•  7000 

1.0709 

-.9505 

8 

•  8000 

1*1069 

-1,0122 

9 

•  9000 

l.lfOl 

-1.0637 

10 

1.0000 

1.1400 

-1.1071 

THt  oXsTawce  F*oh 

ThE  Free  SURFACE  » 

•1*45630  FEET 

*« 

p/RHO.u.A 

PHASE  AHbLE 

1 

•  1000 

.4685 

-.1974 

2 

•2000 

•  2273 

-.3805 

J 

•  3000 

•  1181 

-.3404 

4 

•  4000 

•  0671 

—8747“ 

S 

.5000 

•  Of  §9 

-•7854 

6 

•6000 

.0147 

-.0761 

7 

•  7000 

•  0074 

-#'*85 

8 

•  8000 

•  0087 

-1.0122 

9 

•  9000 

•  0019 

-l*0o37 

10 

1.0000 

.0009 

-1.1071 

RUONTCft  .T  I*  rttv»Ut\CAcS  *  •,>UUtL  ft  CONE  BOiTQM  t*  NU  EXH.) 
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